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Abstract 



We obtain reliable irn scattering amplitudes consistent with experimental 
data, both at low and high energies, and fulfilling appropriate analyticity proper- 
ties. We do this by first fitting experimental low energy {s^^^ < 1.42 GeV) phase 
shifts and inelasticities with expressions that incorporate analyticity and unitar- 
ity. In particular, for the S wave with isospin 0, we discuss in detail several sets of 
experimental data. This provides low energy partial wave amplitudes that sum- 
marize the known experimental information. Then, we impose Regge behaviour 
as follows from factorization and experimental data for the imaginary parts of 
the scattering amplitudes at higher energy, and check fulfillment of dispersion 
relations up to 0.925 GeV. This allows us to improve our fits. The ensuing tttv 
scattering amplitudes are then shown to verify dispersion relations up to 1.42 GeV, 
as well as s — t — li crossing sum rules and other consistency conditions. The 
improved parametrizations therefore provide a reliable representation of pion-pion 
amplitudes with which one can test chiral perturbation theory calculations, pio- 
nium decays, or use as input for CP-violating K decays. In this respect, we find 
[«o°' - ao^? = (0-077 ± 0.008) M"^ and S^''\m%) - s'o\mji) = 52.9 ± 1.6°. 
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1. Introduction 

A precise and unbiased knowledge of the tttt scattering amplitude has become increasingly important in the 
last years. This is so, in particular, because tttt scattering is one of the few places where one has more 
observables than unknown constants in a chiral perturbation theory (ch.p.t.) analysis to one loop, so it 
provides a test of ch.p.t. in this approximation, as well as a window to higher order terms. Beside this, an 
accurate determination of the S wave scattering lengths, and of the phase shifts, provide essential information 
for three subjects under intensive experimental investigation at present, viz., light scalar spectroscopy, pionic 
atom decays and CP violation in the kaonic system. 

In two recent papers, Ananthanarayan, Colangclo, Gasscr and Leutwylerl^^ and Colangelo, Gasser 
and Leutwylerl^l have used experimental information, analyticity and unitarity (in the form of the Roy 
equations!'^!) and, in the second paper, also chiral perturbation theory, to construct a precise tttt scattering 
amplitude at low energy, s^/^ < 0.8 GeV. Unfortunately, however, the analysis of refs. 1, 2 presents some 
weak points. First of all, the input scattering amplitude at high energy (s^/^ > 1.5 GeV) which these authors 
use is incompatible with experimental data^'*! on tttt cross sections and also contradicts known properties of 
standard Regge theory; see a detailed analysis in refs. 5, 6, 7 and Appendix B here. Moreover, the errors these 
authors take for some of their experimental input data are optimistic, as shown in ref. 7. As a consequence 
of all this the tttt amplitudes of ref. 2 do not satisfy well a number of consistency tests, as we show in Sect. 7 
here (see also refs. 6, 7 for more details). 

Some of the shortcomings of the articles in refs. 1, 2, notably incorrect Regge behaviour, are also 
reproduced in the papers of Descotes et al., and Kamihski, Lesniak and Loiseau,'^' who also base their 
analysis in the Roy criuations but, since the errors given by these authors are substantially larger than those 
of ref. 2, their effects are now less pronounced. Therefore, we still need to find reliable pion-pion scattering 
amplitudes compatible with physical data both at high and low energy, as well as to verify to what extent 
such amplitudes agree with ch.p.t. 

In the present paper we address ourselves to the first question; that is to say, we try to find what 
experiment implies for the tttt amplitudes. To avoid biases, we will start by performing fits to experimental 
data on phase shifts and inelasticities, incorporating only the highly safe requirements of analyticity and 
unitarity, in the low energy region s^/^ < 1.42 GeV. In particular, for the SO wave below 0.95 GeV, where the 
experimental situation is confused, we perform a global fit, as well as individual fits to various sets of data. 
These fits are described in Sect. 2 and, in the following Sections, we investigate to which extent the ensuing 
scattering amplitudes are consistent, in particular with high energy information. To do so, we assume Regge 
behaviour, as given in ref. 5 (slightly improved; see Appendix B), above 1.42 GeV: using this we check in 
Sect. 3 fulfillment of forward dispersion relations, for the three independent tttt scattering amplitudes. This, 
in particular, permits us to discriminate among the various sets of phase shifts for the SO wave, leaving only 
a few solutions which are consistent with dispersion relations (and, as it turns out, very similar one to the 
other, as discussed in Sect. 4). 

When dealing with differc^nt data sets one has to weight not only the data on a single experiment 
but one has to take into account the reliability of the experiments themselves. So we have done for many 
waves, where some clearly faulty experimental data have only been considered only to conservatively enlarge 
the uncertainties. Concerning the most controversial SO wave, we have used the very reliable data coming 
from Ki4 and K ^ nn decays; to this we add the results from other experimental analyses of tttt scattering 
available in the literature. We then use forward dispersion to test consistency of the several sets of data. 

The present study should therefore be considered, in particular, as a guideline to the consistency 
(especially with forward dispersion relations) of the various data sets. 

Next, we use these dispersion relations to improve the central values of the parameters of the fits 
given in Sect. 2. The result of such analysis (Sect. 4) is that one can get a precise description for all waves, 
consistent with forward dispersion relations up to s^/^ ~ 0.95 GeV and a bit less so ( < 1.5 cr level) in the 
whole energy range, 2M^ < s^/^ < 1.42 GeV, and even below threshold, down to s^/^ = a/2M^. The greater 
uncertainties affect the SO wave for s^/^ > 0.95 GeV, a not unexpected feature, and, to a lesser extent, the 
P wave above 1.15 GeV. 

In Sect. 5 we verify that the scattering amplitudes we have obtained, which were shown to satisfy 
s — u crossing (by checking the dispersion relations), also verify s — t crossing, in that they satisfy two typical 



- 1 - 



-J. R. PELAEZ AND F. J. YNDURAIN- 



crossing sum rules. In Sect. 6 we use the scattering amplitudes we have determined and the method of the 

Froissart Gribov representation to calculate a number of low energy parameters for P, D and some higher 
waves which, in particular, provides further consistency tests. We also evaluate, in Sect. 7, the important 
quantities [a^'^ — a^^]'^ and S^\mj() — sj^\mj^) for which we find 

- 4'^]' = (0-077 ± 0.008) - 4'' = 52.9 ± 1.6°. 

Also in Sect. 7 we compare our results with those obtained by other authors using Roy equations and ch.p.t. 
However, in the present paper we will not address ourselves to the question of the chiral perturbation theory 
analysis of our tttt amplitudes. 

Our paper is finished in Sect. 8 with a brief Summary, as well as with a few Appendices. In 
Appendix A, we collect the formulas obtained with our best fits. In Appendix B we give a brief discussion 
of the Reggc formulas used; in particular, we present an improved evaluation of the parameters for rho 
exchange. Appendix C is devoted to a discussion of the shortcomings of experimental phase shift analyses 
above ^ 1.4 GeV, which justifies our preference for using Regge formulas in this energy region. 

We end this Introduction with a few words on notation and normalization conventions. We will here 
denote amplitudes with a fixed value of isospin, say /, in channel s, simply by F'^^\ //^''i we will specify the 
channel, F^-^'\ when there is danger of confusion. For amplitudes with fixed isospin in channel t, we write 
explicitly F^^*\ 

For scattering amplitudes with well defined isospin in channel s. Is, we write 
F'^^'\s,t) =2x J2 (2« + l)Pi(cos0)//^=^(s), /, = even, 

i=even 

F^''\s,t) = 2 X ^{21 + l)Pi{cose)fl'^\s), I, = odd, 

i=odd (1.1a) 

//''(^)=^r, r=sin^F^(s)e-l^'(^). 



The last formula is only valid when only the clastic channel is open. When inelastic channels open this 
equation is no more valid, but one can still write 



Ms) = 



1 



2i 



(1.1b) 



The factor 2 in the first formulas in (1.1a) is due to Bose statistics. Because of this, even waves only exist 

with isospin / = 0, 2 and odd waves must necessarily have isospin / = 1. For this reason, we will often omit 
the isospin index for odd waves, writing e.g. /i, /a instead of f\^\ Another convenient simplification 

that we use here is to denote the tttt partial waves by SO, S2, P, DO, D2, F, etc., in self-explanatory notation. 

The quantity rji, called the inelasticity parameter for wave Z, is positive and smaller than or equal 
to unity. The elastic and inelastic cross sections, for a given wave, are given in terms of 5i and rii by 



el. 



}• 



CT; ■, a^' are defined so that, for collision of particles A, B (assumed distinguishable). 



4^2 2si/2 
crtot. = TYT^T 7 —j- 22(21 + 1) [af- + a, 



inel. 



(1.2) 



(1.3) 



When neglecting isospin violations (which we do unless explicitly stated otherwise) we will take the 
convention of approximating the pion mass by = JTiTri — 139.57 MeV. We also define scattering lengths, 

by 



a^^^ , and effective range parameters, 



TT 



(1.4) 



-2- 



-THE PION-PION SCATTERING AMPLITUDE- 

2. The TTTT scattering amplitudes at low energy (s^/^ < 1.42 GeV) 

Wc present in this section parametrizations of the S, P, D, F and G waves in tttt scattering obtained fitting 
experimental data at low energy, s^/^ < 1.42 GeV, which will provide a representation of the experimental 
TTTT scattering amplitudes in this energy range: this may be considered as an energy-dependent phase shift 
analysis. Above 1.42 GeV, wc will use the Rcggc expressions obtained in ref. 5, which wc reproduce and 
discuss in Appendix B in the present paper. In following Sections we will verify to which extent the scattering 
amplitudes that one finds in this way are consistent with dispersion relations or crossing symmetry, within 
the errors given. 

Before entering into the actual fits, a few words are due on the choice of the energy 1.42 GeV as 
limit between the regions where we use phase shift analyses and Regge representations. Experimental phase 
shift analyses exist up to ^ 2 GcV. However, as is known, phase shift analyses become ambiguous as soon 
as inelastic processes become important. As we will discuss in some detail in Appendix C, the existing 
experimental phase shift analyses become suspect for energies above ~ 1.4 GeV: in particular, they produce 
cross sections that deviate from experimentally measured total cross sections. 

As for the higher energy range, we have shown in ref. 5 that a Regge description of the imaginary 
part of the tttt scattering amplitudes agrees well with tttt cross section data (and, through factorization, also 
with ttN and NN data) for kinetic energies above ^ 1 GeV; for tttt scattering, down to s^^'^ ~ 1.4, or even 
to 1.3 GeV, depending on the process. We have chosen the limiting energy to be 1.42 GeV as a reasonable 
compromise; but one could have taken lower limits, or slightly higher ones, as well. In fact, it is not possible 
to go below 1.4 GeV with a Regge description for the It = I amplitude or for tt^tt" scattering, since both 
contain the isospin zero amplitudes, which show two resonances just below that region [the /2(1270) and the 
/o(1370)]. But it is possible to choose a lower junction point between the phase shift analyses and the Regge 
formulas for 7r°7r+ scattering. The influence of this is negligible for tt^tt^, as wc will show in Subsect. 3.1.2, 
provided one stays in the range 1.32 GeV < s^^^ < 1.42 GeV, with the larger values favoured. 

We now turn to the parametrizations. In writing them we will use the requirements of analyticity 
and elastic unitarity. Extra information has been added to help stabilize the fits in those channels where the 
low energy data is scarce; in particular, information on Adler zeros or scattering lengths. For the latter we 
impose values obtained from their Froissart-Gribov representation, but with very large error bars to cover the 
difiierent values available in the literature. The effect on this of the high energy representation is negligible. 

The method used to take into account unitarity and analyticity will be the standard one of the 
eff'ective range formalism, supplemented by a conformal expansion. To be precise, for a given partial wave 
fl^\s) we write, for any value (complex or real) of s, 

/n 2sl/2 1 1.2/ 1 

f('\s) = — = = - = (2 1) 

' nk 2sV2fc-2i-l^W(s) -i TT ^W(s)_ifc2m/2sl/2" 

The effective range function ^i^\s) is real on the segment < s < sq, but it will be complex above sq, and 
also for s < 0. Here sq is the energy squared above which inelastic processes are nonnegligiblc. Using only 
the requirements of causality and conservation of probability, it can be shown that <?[^''(s) is analytic in the 
complex s-planc, cut from — oo to and from Sq to +oo; see Fig. 1. 

We can profit from these analytical properties as follows. We map the cut plane into the unit disk 
(as in Fig. 1) by means of the conformal transformation 



w{s) 



y/s + ^ySQ - S 



The properties of analyticity and elastic unitarity of fj;^^ (s) are then strictly equivalent to uniform conver- 
gence (in the variable w) in the disk \w\ < 1 of the series 

oo 
n=0 
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FlGURE 1. 

The mapping s w. 



w 



s 



This series thus converges, in the variable s, uniformly in the interior of the whole cut plane; in general very 
quickly. It will turn out that we will need only two or three terms in the expansions for the partial wave 
amplitudes. 

Before starting with the actual fits, we here say a few words on the values of so (i-e., the energy at 
which we consider inelasticity not negligible) that we will take for the various waves. For the SO wave, the 
KK channel is strongly coupled, so here sq = 4m^. For other waves the KK channel is weakly coupled. In 
the cases where we have sufficiently precise data (that is, for the P wave and the low energy S2 wave), we 

1/2 

will take Sq — 1.05 GeV which is, approximately, the 27rp threshold. 

For the DO wave, and for the S2 wave at intermediate energy, 1.05 GeV < 5^/2 < 1.4 GeV, inelas- 
ticity is detected. It is small (but not negligible^) while, unfortunately, the data arc not precise enough to 
perform a fully consistent analysis. We then follow a different strategy: we fit the experimental phase shifts 
using formulas that neglect inelasticity below 1.45 GeV (which is, approximately, the 2p threshold). The 
experimental inelasticity between 1 GeV and 1.42 GeV is then added by hand. 

Finally, for the D2, F waves, for which experiments detect no inelasticity below 1.42 GeV, Sq is taken 
equal to 1.45 GeV. 

2.1. The P wave 

2.1.1. Parametrization of the P wave below 1 GeV 

We will consider first the P wave for tttt scattering below 1 GeV, to exemplify the methods, because it is the 
one for which more precise results are obtained. We start thus considering the region of energies where the 
inelasticity is below the 2% level; say, sj^^ < 1.05 GeV . 

We then expand ^i(s) in powers of w, and, recxpressing w in terms of s, the expansion will be 
convergent over all the cut s-plane. Actually, and because we know that the P wave resonates at s = Mp, 
it is more convenient to expand not ^i(s) itself, but V(s) given by 



In terms of ^i{s) we thus find the expression for the cotangent of the phase shift, keeping two terms 
in the expansion. 



#i(s) = (s-M2)V.(s)/4; 



(2.2a) 



so we write 



i){s) = {Bo + Biw + ■ ■ ■} ; 



(2.2b) 




Sq/^ = 1.05 GeV. 



(2.3) 
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Mp, Bo, Bi are free parameters to be fitted to experiment. In terms of tp we have, for the rho width, 
and the scattering length, ai, 



3 



^ -1 ^ 1 ^^-^^ 

""^ ~ 4M„#i(4M2) ~ M^(M2 - 4M2)V(4M2) " 

The values Bq = const., Bi>i = would correspond to a perfect Breit-Wigner. Actually, it is known 

that the p deviates from a pure Brcit- Wigner and for a precision parametrization two terms, Bq and Bi, 
have to be kept in (2.3). Note that the parametrization holds not only on the physical region 4M^ < s < sq, 
but on the unphysical region < s < 4M^ and also over the whole region of the complex s plane with 
Im s ^ 0. The parametrization given now is the one that has less biases, in the sense that no model has 
been used: we have imposed only the highly safe requirements of analyticity and unitarity, which follow from 
causality and conservation of probability. 

The best values for our parameters arc actually obtained not from fits to tttt scattering data, but from 
fits to the pion form factor. 1^1 This is obtained from eir scattering, from e^e~ tt+tt" and from r vittt 
decay, where we have a large number of precise data, and pions are on their mass shells. Including systematic 
experimental errors in the fits, and including in the fit also the constraint ai = (38 ± 3) x 10~^ foi' the 
scattering length, we have 

fromTT+TT": Bq = 1-074 ± 0.006, Bi = 0.13 ± 0.04; Mp = 774.0 ± 0.4 MeV , . 

from7r+7r°: Bq = 1-064 ± 0.006, Si = 0.13 ± 0.04; Mp = 773.1 ± 0.6 MeV . ^ ' 

The fit is excellent; one has /d.o.f . = 245/(244 — 13), and the slight excess over unity is due to the 
well-known incompatibility of OPAL and other data for r decay at low invariant mass. 

The corresponding values for the rho width, and P wave scattering length and effective range pa- 
rameter, are 

ai = (37.8 ± 0.8) x IQ-^M'^, 6i = (4.74 ± 0.09) x IQ-^M'^, Fpo = 146.0 ± 0.8; from tt+tt" 

ai = (37.8 ± 0.8) x lO'^M-^, h = (4.78 ± 0.09) x IQ-^M'^, = 147.7 ± 0.7; from 7r+7r°. ^'^'^^^ 

Although the values of the experimental tttt phase shifts were not included in the fit, the phase shifts that 
(2.5a) implies are in very good agreement with them, as shown in Fig. 2. 

If neglecting violations of isospin invariance, one would expect numbers equal to the average of both 
sets. If wc also increase the error, so as to take into account that due to isospin breaking, by adding linearly 
half the difference between both determinations (2.5a), we find 

Bo =1-069 ±0.011, Bi =0.13 ±0.05, Mp = 773.6 ± 0.9, 

ai = (37.6 ± 1.1) X lO-^M-^, hi = (4.73 ± 0.26) x lO'^M"^; ^'^'^^ 
the corresponding parametrization we take to be valid up to s^^^ = 1 GeV. 



2.1.2. The P wave for 1 GeV < s^'^ < 1.42 GeV 

In the range 1 GeV < s^/^ < 1.42 GeV one is sufficiently far away from thresholds to neglect their influence 

(the coupling to KK is negligible). A purely empirical parametrization that agrees with the data in Pro- 
topopescu et al.[^°l and Hyams et al.,^^^"' up to 1.42 GeV, within errors, is obtained with a linear fit to the 
phase and inelasticity: 

^i(s) =Ao + Ai(vW^- 1), ?7i(s) = 1 - e(vW^- 1); 
e= 0.30 ±0.15, Ac = 2.69 ±0.01, Ai = l.l±0.2. 

Here s = 1 GeV^. The value of Aq ensures the agreement of the phase with the value given in the previous 
Subsection at s = s = 1 GeV^. This fit is good (Fig. 2). 
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Figure 2. The phase shifts of solution 1 from Protopopescu et al.,'^"' Hyams et al.'^^"' and Es- 
tabrooks and Martin'^^"' compared with the prediction with the parameters (2.6) (solid line below 
1 GeV). We emphasize that this is not a fit to these experimental data, but is obtained from the pion 
form factor.'^' The error here is like the thickness of the line. 
Above 1 GeV, the dotted line and error (PY) are as follows from the fit (2.7). 

It should be remarked, however, that there are other solutions for the P wave above 1.15 GeV. In 
the first analysis of Hyams et al.,[^^°l that we have followed here, a resonance occurs with a mass ~ 1.6 GeV, 
and its effect is only felt above ~ 1.4 GeV; but in the second analysis by the same group, Hyams et al.,[^^^l a 
broad, highly inelastic resonance (actually, more a spike than a resonance) with a mass ~ 1.35 GeV appears, 
or does not appear, depending on the solution chosen. Finally, the Particle Data Group (based mostly on 
evidence from e"'"e~ annihilation and r decay) report a resonance at 1.45 GeV: one has to admit that the 
situation for the P wave above 1.15 GeV is not clear. We will return to this later. 

2.2. The S waves 

2.2.1. Parametrization of the S wave for 1 = 2 

We consider three sets of experimental data. The first corresponds to solution A in the paper by Hoogland et 
al.,[^^l who use the reaction 7r+p ■K^n'^n. We will not consider the so-called solution B in this paper; while 
it produces results similar to the other, its errors are clearly underestimated. The second set corresponds 
to the work of Losty et al.,'"'^^] who analyze instead n^p — > tt^tt^A. The third set are the data of Cohen 
et al.,W who consider Tr~n — > TT~i:~p. These three sets represent a substantial improvement over other 
determinations; since they produce two like charge pions, only isospin 2 contributes, and one gets rid of the 
large SO wave contamination. However, they still present the problem that one does not have scattering of 
real pions. 

For isospin 2, there is no low energy resonance, but /g (s) presents the feature that a zero is expected 
(and, indeed, confirmed by the fits) in the region < s < 4Af^. This zero of f^\s) is related to the so-called 
Adler zeros and, to lowest order in chiral perturbation theory, occurs at s = 2z| with Z2 = M^. We note 
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FlGURE 3. Continuous line: The 1 = 2, S-wa.ve phase shifts, and error bands, corresponding to 
(2.7), (2.8) and (2.9), denoted by PY. Dotted line: the fit from ref. 6. The experimental points are 
from Losty et al. (black squares), Hoogland et al.,'^^' solution A (black dots), and Cohen et al.W (black 
triangles). We have also included data from Durusoy et al. and from Solution B of Hoogland et al., 
although they were not used in the fits. The dashed line, which lies below our fit, is the S2 phase of 
Colangelo, Gasser and Leutwyler,''^' with error band attached. 



that, unlike the corresponding zero for the SO wave, 2z| is inside the region where the conformal expansion 
is expected to converge. 

In ref. 6, we neglected inelasticity below 1.45 GeV, and fitted all experimental data, for s^/^ < 
1380 MeV. A more precise determination (and, above all, with more realistic error estimates) is obtained 
if we realize that inelasticity is detected experimentally above the 2ttp threshold, s^/^ — 1.05 GeV, so we 
should fit separately the low and high energy regions (Fig. 3). 

In the low energy region, we fix Z2 = and fit only the low energy data, s^^^ < 1.0 GeV; later, in 
Sect. 4, we will allow Z2 to vary. We write 

cot4^)(g)^^^^ (i3o + i?i ^'^^ j, Z2^M.; 1.05 GeV. (2.8a) 

ZIt s — I y's + V s — s J 

Then we get /d.o.f. = 13.0/(25 - 2) and 

Bo = -80.4 ±2.8, Bi -73.6 ± 12.6; = (-0.052 ± 0.012) Af^^ 6^^^ = (-0.085 ± 0.011) Af-^. 

(2.8b) 

For the high energy region we neglect the inelasticity below 1.45 GeV, and then add inelasticity by 
hand. We consider two extreme possibilities: fitting the whole range, or fitting only high energy data (s^/^ > 
0.91 GeV), requiring agreement of the central value with the low energy determination at s^/^ = 1 GeV. We 
accept as the best result an average of the two, and thus have 



.a(2)^ A g'^' , p \/^-\/^ 



2k S-2A/2 t V^+Vi^7^J' (2.9) 

^1.45 GcV; Bo = -123 ±6, Bi = -118±14, 
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and we have enlarged the errors to cover both extreme cases. We will not consider varying the position of 
the Adler zero for this high energy piece. 

The inelasticity may be described by the empirical fit 

„(2)/,x_/l-e(l-s/s)3/2, e = 0.18 ±0.12; s > s = (1.05 GeV)^; 

'^o ~ 1 1, s<s = (1.05 GeV)2. ^^-^"^ 

These formulas are expected to hold from s^/^ = 1.0 GeV up to 1.42 GeV. As shown in Fig. 3, both 
determinations, (2.8-9) and that in ref. 6 are almost overlapping (their error bands actually overlap). 

2.2.2. Parametrization of the S wave for 7 = below 0.95 GeV (global fit) 

The S wave with isospin zero is difficult to deal with. Here we have a very broad enhancement, variously 
denoted as e, a, /o, around s^^^ = Mo ~ 800 MeV. We will not discuss here whether this enhancement is 
a bona Bdes resonance; we merely remark that in all experimental phase shift analyses Sq^\s) crosses 90° 
somewhere between 600 and 900 MeV; we define /kq as the energy at which the phase equals 90°. Moreover, 
we expect a zero of f^°\s) (Adler zero), hence a pole of the effective range function <?q°'(s), for s = ^Zq 
with ^Zq in the region < s < 4M^. Chiral perturbation theory suggests that zq — M„. 

1/2 ~ 

We can distinguish two energy regions: below Sq = 2mK we are under the KK threshold. Between 

sj^^ and s^/"^ ^ 1.2 there is a nonegligiblc coupling between the KK and tttt channels and the analysis 
becomes very unstable, because there is little information on the process tttt — » KK and even less on 
KK — > KK. We will present later an empirical fit in the region of energies around and above 1 GeV, and 
we will now concentrate in the low energy region. 

For the theoretical formulas we impose the Adler zero at s = ^Af^ (no attempt is made to vary this 

for the moment; see Sect. 4), and a zero of cot(5g'''(s) at s = /Xg, /kq a free parameter. Then we map the s 
plane, cut along the left hand cut (s < 0) and the KK cut, writing 

„l/2 j^-2 ,,2 „ , ^ 

coiCis) = ^{Bo + B,w{s) + ...}, 

s — ^Zq Mo ^ ^ 

and 

, \/s — \/So — S ^ 2 

^(^) = 1^, — / So = 4m^ 

y/S + VSo - S 

(we have taken niK = 0.496 GeV). 

This parametrization does not represent fully the coupling of the KK channel and we will thus only 
take it to be valid up to s^/^ = 0.95 GeV. 

On the experimental side the situation is still a bit confused, although it has cleared up substantially 
in the last years thanks to the experimental information on K14 and K2Tr decays. The information we have on 
this SO wave is of three kinds: from phase shift analysis in collisionsl^^'^^"' irp ttttN, A; from the decayl^^l 
K14: and from K2n decays. The last gives the value of the combination sj^^ — Jg^'' at s^^"^ = tuk] the decay 
Ki4^ gives (5g"'' — 5i at low energies, s^/^ < 380 MeV. If using recent K2.n information, I^^l combined with 
older determinations, and with the 1 = 2 phase obtained in the previous subsection, one finds the phase 

Sl"\m]^)=43.3±3° . (2.11) 

We will here include in the fit the low energy data from K14 decay, ^ shown in Fig. 4, and we impose 
the value of ^g°^(?7i|-) from K2Tr given in (2.11). The main virtue of these K decay data is that they refer 

to pions on their mass shell; but, unfortunately, this leaves too much room at the upper energy range, 
s^/^ > 0.6 GeV. If we fit only K decay data we can only use one parameter Bq in the conformal expansion: 



As a technical point, wo mention tliat we have increased by 50% the error in the point at highest energy, s^^^ = 
381.4 MeV, from the Ke4 compilation of Pislak et al., because this experimental value represents an average over 
a long energy range that extends to the edge of phase space. 
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FlGURE 4. The / = 0, 5- wave phase shifts and error bands corresponding to Eq. (2.14) (PY, contin- 
uous line). Also shown (black dots) are the points from Ku and K^-n decays, and only the high energy 
data included in the fits, as given in Eq. (2.13). The dashed line is the solution of Colangelo et al.'^' 



if including another parameter, spureous minima would appear. We get a good fit, albeit with rather large 
errors: 

2^0 /^o 

y2 5.7 (2.12) 

So =18.5 ±1.7, = 766 ± 95 MeV; = ——. ^ ' 

d.o.i. 12 — 2 

4"^ = (0.22 ± 0.02) X Af-i; 5^^\mK) = 43 ± 5° . 

To improve on this we have to add further data (and one more parameter Bi). To do so, we can 
follow two different procedures. We can add to the K^i, K2-n data various sets of experimental phase shifts, 
fitting each set individually; this we will do in Subsect. 2.2.3. Or we can follow what we consider to be the 
best procedure: we combine in the fit data from various experiments at energies above 0.8 GeV, which we 
will call a global fit. The reason to choose only data above 0.8 GeV is that, in the region between 0.81 GeV 
and 0.97 GeV, the more relevant experimental results have overlapping error bars, something that does not 
happen at other energies (see Fig. 5). By combining several sets we may expect to average out systematic 
errors, at least to some extent. 

At high energy we thus include the following sets of data: first of all, the values 

6^°^ (0.870^ GeV^) = 91 ± 9° ; (0.910^ GeV^) = 99 ± 6° ; 

(2.13a) 

6^°> {0.935^ GeV^) =109 ± 8° ; S^°\o.965'^ GeV^) = 134 ±14°. 

These points are taken from solution 1 of Protopopescu et al.'^*^' (both with and without modified moments), 
with the error increased by the difference between this and Solution 3 data in the same reference. We will 
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FlGURE 5. The 1 = 0, S-wave phase shifts and error band in the whole energy range as given by our 
fits (PY), Eqs. (2.f4), (2.f5). The experimental points from K14 and K2tt decays are not shown. 

also mclude in the fit the data, at similar energies, of Grayer et al.:!^^"' 

4''^(0.9122 GeV^) =103±8°; 5(o)(o.9292 GeV^) = 112.5 ± 13° ; 

S^°\0.952'^ GeV^) = 126 ± 16° ; 6^°\o.970'^ GeV^) ^ 141 ± 18° . 

The central values are obtained averaging the solutions given by Grayer et al., except solution E, and the 
error is calculated adding quadratically the statistical error of the highest point, the statistical error of the 
lowest point (for each energy) and the difference between the central value and the farthest point. 

Finally, we add three points between 0.8 and 0.9 GeV obtained averaging the s-channel solution of 
Estabrooks and Martin and solution 1 of Protopopescu et al., which represent two extremes. The error is 
obtained adding the difference between these two in quadrature to the largest statistical error. In this way 
we obtain the numbers, 

4"^(0.8102 GeV^) =88±6°; ^"^(O.SSO^ GeV^) = 92 ± 7° ; sj°\o.850^ GeV^) ^ 9i± 6° . (2.13c) 

In spite of the generous errors taken, it should be noted that these data could still contain systematic 
errors, beyond those taken into account in (2.13), which may contaminate the results of the fit at the higher 
energies. In Sect. 4 we will use dispersion relations to improve the parametrization of this SO wave. 

We now have enough constraints to include two parameters Bi in the expansion of cot ^q""* . We find 



Bo =21.04, Bi = 6.62, /io = 782 ± 24 MeV; xVd-O.f. = 15.7/(19 - 3). 
4°^ =(0.230±0.010)M-\ = (0.268 ±0.011)1/^3; ^("'(m^^) 41.0° ±2.1° 



s 



(2.14a) 
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this fit (shown in Fig. 4) we take to be valid for s^/^ < 0.95 GeV. The errors of the Bi are strongly correlated; 
uncorrected errors are obtained if replacing the -B, by the parameters x, y with 

Bo=y-x; Bi = 6.62 - 2.59a;; ?/ = 21.04 ± 0.70, a; = 0±2.6. (2.14b) 
2.2.3. Parametrization of the S wave for 7 = below 0.95 GeV (individual fits) 

In this Subsection we summarize, in Table 1, the results of fits to data from K14 and K2Tr decays including 
also, individually, data from various sets of phase shift analyses. The method of fit is like that used in the 
previous Subsection; in particular, we have fixed the Adler zero at zq = in all these fits. 





Bo 




/io (MeV) 




PY, Eq. (2.14) 


21.04 


6.62 f''' 


782 ± 24 


0.230 ± 0.010 


K decay only 


18.5 ± 1.7 


= 


766 ± 95 


0.218 ±0.021 


K decay data 
+ Grayer, B 


22.7 ± 1.6 


12.3 ±3.7 


858 ± 15 


0.246 ± 0.042 


K decay data 
+ Grayer, C 


16.8 ±0.85 


-0.34 ±2.34 


787 ±9 


0.236 ± 0.023 


K decay data 
+ Grayer, E 


21.5 ±3.6 


12.5 ± 7.6 


1084 ± 110 


0.26 ± 0.05 


K decay data 
+ Kaminski 


27.5 ±3.0 


21.5 ± 7.4 


789 ± 18 


0.25 ±0.10 


K decay data 
+ Grayer, A 


28.1 ± 1.1 


26.4 ± 2.8 


866 ±6 


0.29 ± 0.04 


K decay data 
+ EM, s— channel 


29.8 ± 1.3 


25.1 ±3.3 


811 ±7 


0.27 ±0.05 


K decay data 
+ EM, t-channel 


29.3 ± 1.4 


26.9 ± 3.4 


829 ±6 


0.27 ±0.05 


K decay data 
+ Protopopcscu VI 


27.0 ± 1.7 


22.0 ±4.1 


855 ± 10 


0.26 ± 0.05 


K decay data 
+ Protopopescu XII 


25.5 ± 1.7 


18.5 ±4.1 


866 ± 14 


0.25 ±0.05 


K decay data 
+ Protopopescu VIII 


27.1 ±2.3 


23.8 ± 5.0 


913 ± 18 


0.27 ±0.07 



PY, Eq. (2.14): our global fit, Subsect. 2.2.2. The rest are fits to either K decay data alone, or combining these with 
various tttt scattering data sets. Grayer A, B, C, E: the solutions in the paper of Grayer et al.'^^"' (solution D only 
concerns data above 1 GeV). EM: the solutions of Estabrooks and Martin. Kaminski: the papers of Kaminski 
et al.'^^*^' Protopopescu VI, XII and VIII: the corresponding solutions in ref. 10. Solutions A, B and C of Grayer 
et al., as well as both solutions EM and Kaminski et al. are from energy-independent analyses. The rest are from 
energy-dependent evaluations. '■^^ We do not give errors here as they are strongly correlated; cf. Eq. (2.14b). 

Table 1 



A few comments arc in order. First of all, wc note that the solution E of Grayer ct al., as well as 
what one finds with only K decay data, have very large errors. Moreover, solution E of Grayer et al. only 
contains eight points which clearly lie below all other experimental data; see Fig. 5. Second, it is seen that 
the values of the parameters in the first five fits (which, as will be shown in Sect. 4, are the more reliable 
ones) cluster around our solution, labeled PY, Eq. (2.14) in Table 1. This is as should be expected, since 
our global solution sums up information from various experimental sets. Third, unlike our global solution, 
which includes systematic errors, the errors given in Table 1 for the other fits are only the statistical ones. 
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Finally, we remark that all fits other than the first six have parameters that differ a lot from those obtained 

fitting only K decay data. Since this last fit is already very good, this makes the lower six fits in Tabic 1 
suspect. Indeed, and as we will see, they lead to scattering amplitudes that do not satisfy well dispersion 
relations. 



2.2.4. The / = S wave between 950 MeV and 1420 MeV 



The description of pion-pion scattering above the KK threshold would imply a full two-channel formalism. 
To determine the three independent components of the effective range matrix <1>22 and <?i2, one 

requires measurement of three cross sections. Failing this, one gets an indeterminate set, which is reflected 
very clearly in the wide variations of the effective range matrix parameters in the energy-dependent fits of 
Protopopescu et al.[^°l and Hyams et al.I^^^l 

The raw data themselves are also incompatible; Protopopescu et al. find a phase shift that fiattens 
above s^^"^ ~ 1.04 GcV, while that of Hyams et al. or Grayer et al. continues to grow. This incompatibility 
is less marked if we choose the solution with modified higher moments by Protopopescu et al. (Table Xlll 
there). The inelasticities are compatible among the various determinations from tttt scattering, including 
that of Hyams et al.,[^^^! and Kaminski ct aL.t^^'^I although the errors of Protopopescu ct al. appear to be 
much underestimated. They are, however, systematically, well above what one finds!"''' from TTTT KK 
scattering: see Fig. 6. 

We will here give a purely empirical fit, using the tttt data. We write 



cotJ^^(.)=co ^, 

,(0) = 1 - (er + e^M) M' = 1.5 GeV (fixed). 



(2.15a) 



The fit to the inelasticity gives 

ei=6.4±0.5, €2 = -16.8 ± 1.6; xV(d-o.f.) = 0.7. (2.15b) 

This result is driven by the data of Protopopescu et al., whose accuracy is, unfortunately, much overestimated: 
see Fig. 6. 

If, instead of fitting 7]^^^ to the tttt data of Protopopescu et al. and Hyams et al. we had fitted the 
data!"''] from TTTT — > KK scattering (shown in Fig. 6), we would have found values for the Cj much smaller 
than what was given in (2.15b): 

ei = 2.4±0.2, 62 = -5.5 ±0.8; xV(d-o.f.) = 1.3. (2.15b') 

We have checked that the influence of using (2.15b) or (2.15b') on the dispersion relations and other 
evaluations of integrals, to be considered later, is minute, for energies below 0.95 GeV. This is because the 
inelasticity affects little the imaginary part of the partial wave (on the average). Above 1 GeV, if we took 
the Ci following from tttt — > KK, Eq. (2.15b'), the dispersion relations would be slightly better fulfilled; see 
Subsect. 4.1, at the end. In spite of this, we stick to (2.15b). Taking rjj^'^ from one set of experiments and 
(5q°' from another (incompatible with the first) would be an inconsistent procedure. 

To fit (5g°^ we also require it to agree with the low energy determination we found in Subsect. 2.2.2 
at = 0.95 GeV. If we include the data of Protopopescu et al. in the fit we find a poor fit with 
/d.o.f. = 39/(14 - 2) and the parameters 

Co = 1.72 ± 0.08, M, = 930 MeV, Mf = 1340 MeV . 

The error in cq, corresponding to Scr, is purely nominal. If we keep Mg fixed and remove the data of 
Protopopescu et al., we get a good x^/d.o.f., and now 

Co = 0.79 ± 0.25, Mf = 1270 MeV . 
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Figure 6. Fit to the 1 = 0, 5- wave inelasticity and phase shift between 950 and 1400 MeV. Data 
from refs. 10, 11. The difference between the determinations of ry^"' from tttt — > tttt (PY from data) and 
from TTTT — > KK (PY alternative) is apparent here. 



If we want to be compatible with the data of Hyams et al., we must increase the errors. We then take the 
numbers 

Co = 1.3 ± 0.5, Mf = 1320 ± 50 MeV, = 930 MeV (fixed). (2.16) 

The fit is shown in Fig. 6. 

We emphasize again that these are purely empirical fits and, moreover, they are fits to data which 
certainly have uncertainties well beyond their nominal errors, as given in given in (2.15b) and (2.16); some- 
thing that is obvious for rj^"' from Fig. 6. It follows that relations such as dispersion relations in which the 
SO wave plays an important role will be unreliable for energies near and, especially, above KK threshold 
(below these energies, however, both (2.15b) and (2.15b') give very similar results). In fact, we will check 
that mismatches occur when s^/^ > 0.95 GeV. A sound description of the SO wave for > 0.95 GeV m 
TTTT scattering would require more refined parametrizations and, above all, use of more information than just 
TTTT experimental data, and lies beyond the scope of the present paper. 

2.3. The D waves 

The D waves cannot be described with the same accuracy as the S, P waves. The data are scanty, and have 
huge errors. That one can get reasonable fits at all is due to the fact that one can use low energy information 
from sum rules; specifically, we will impose in the fits the values of the scattering lengths that follow from 
the Froissart-Gribov representation. Note that this is not circular reasoning, and it only introduces a small 
correlation: the Froissart-Gribov representation for the DO, D2, F waves depends mostly on the SO, S2 and 
P waves, and very little on the DO, D2, F waves themselves. 
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2.3.1. Parametrization of the 7 = 2 D wave 



For isospin equal 2 we would only expect important inelasticity when the channel tttt pp opens up, so we 
will take the value sq = 1.45^ GeV^ ~ AMp for the energy at which elasticity is not negligible. 

But life is complicated: a pole term is necessary to get an acceptable fit down to low energy since wc 
expect S^^ to change sign near threshold. The experimental measurements (Losty et al.; Hoogland et al.I^^l) 
give negative and small values for the phase above some 500 MeV, while chiral perturbation calculations and 
the Froissart Gribov representation indicate a positive sc;attcring length. If we want a parametrization that 
applies down to threshold, we must incorporate this zero of the phase shift. What is more, the clear inflection 
seen in data around 1 GeV implies that we will have to go to third order in the conformal expansion. So we 
write 

cot <5f (a) = ^ {^0 + B,w{s) + B^wisf} ^^^ff^,^ _ ^ (2.17a) 
with A a free parameter fixing the zero and 



x^-x^^ ,1/2 = 1450 MeV. 

Since the data we have on this wave are not accurate (cf. Fig. 7) we have to include extra information. To 
be precise, we include in the fit the value of the scattering length that follows from the Froissart-Gribov 
representation, 

= (2.78 ±0.37) X 10-^M-^ 
but not that of the efi^ective range parameter, 

6(2) = (-3.89 ±0.28) X 10"^ M^"^ 

(see below. Sect. 6). 

We get a mediocre fit, /d.o.f. = 71/(25 — 3), and the values of the parameters are 
Bo = (2.4±0.3) X 10^ = (7.8±G.8) X 10^, B2 = (23.7 ± 3.8) x 10^, Zi = 196 ± 20 MeV . (2.17b) 

We have rescaled the errors by the square root of the x^/d.o.f. 

The fit, which may be found in Fig. 7, returns reasonable numbers for the scattering length and for 

(2) 

the effective range parameter, 63 '■ 

af^ = (2.5 ± 0.9) X 10"^ M^'^; bf^ = (-2.7 ± 0.8) x 10"^ M^"^. (2.18) 

Although the twist of S^\s) at s^/^ ~ 1.05 GeV is probably connected to the opening of the 2'!rp 
channel, wc neglect the inelasticity of the D2 wave, since it is not detected experimentally. This, together 
with the incompatibility of the various sets of experimental data and the poor convergence of the conformal 
series, indicates that the solution given by (2.17) is, very likely, somewhat displaced with respect to the 
"true" D2 wave at the higher energy range (say for s^/^ > 0.7 GeV). In fact, the values of the parameters 
will be improved in Sect. 4 with the help of dispersion relations; the D2 phase shift one finds by so doing is 
slightly displaced with respect to that following from (2.17), as shown in Fig. 7. 
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FlGURE 7. Fits to the 1 = 2, D-wave phase shift. Continuous Hue (PY) with Eq. (2.18). Dashed hue: 
PY, after improving with dispersion relations (Sect. 4). Dotted hue: the fit used in refs. 1, 2 (ACGL, 
CGL). Also shown are the data points of Losty et al., from solution A of Hoogland et al.,'^^' from Co- 
hen et al.'*' and Durusoy et al.,'^^' the last nont included in the fit. 



2.3.2. Parametrization of the / = D wave 



The D wave with isospin in tttt scattering presents two resonances below 1.7 GeV: the /2(1270) and 
the /2(1525), that we wiU denote respectively by /2, f2- ExperimentaUy, Ff.^ = 185 ± 4 GeV and //^ — 
76 ± 10 GeV . The first, /a, couples mostly to tttt, with small couphngs to KK (4.6 ± 0.5 %), 47r (10 ± 3 %) 
and 7777. The second couples mostly to KK, with a small coupling to 7777 and 27r, respectively 10 ± 3 % and 
0.8 ±0.2%. This means that the channels tttt and KK are essentially decoupled and, to a 15% accuracy, we 
may neglect inelasticity up to s ~ 1.45^ GeV^. 

There are not many experimental data on the D wave which, at accessible energies, is small. So, 
the compilation of 6^^^ phase shifts of Protopopescu et al.'^^l gives significant numbers for 5'"2^ only in the 
range 810 MeV < s^/"^ < 1150 MeV. In view of this, it is impossible to get accurately the D wave scattering 
lengths, or indeed any other low energy parameter, from this information: so, we will include information 
on 02"'' to help stabilize the fits. 

We take the data of Protopopescu et al.'^^l and consider the so-called "solution 1", with the two 
possibilities given in Table VI and Table XIII (with modified higher moments), in the range mentioned before, 
gi/^ = 0.810 GeV to 1.150 GeV. The problem with these data points is that they are certainly biased, as 
indeed they are quite incompatible with those of other experiments. We can stabilize the fits by fitting also 
the points^ of Estabrooks and Martin, and imposing the value of the width of the /a resonance, with the 
condition Ff^ — 185 ± 10 MeV, as well as the value of the scattering length. We write 

1/2 

cotd!f \s) = 1^ {Ml - s) AU^ [Bo + Biwis)} (2.19a) 
^ For these data we arbitrarily take a common error of 10%. 
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FlGURE 8. Fit to the / = 0, _D-wave phase shift. Also shown are the data points from solution 1 of 
Protopopescu ct al. (open circles) and data of Estabrooks and Martin (black dots). The dotted lines 
mark the /2(1270) resonance, whose location and width were included in the fit. 



and 

w(s) = s^^ = 1450 MeV; M/, = 1275.4 MeV . 

V s + vso - s 

We find 

5o = 23.6 ±0.7, Si =24.7 ±1.0. (2.19b) 

The fit is not good in that wo get /d.o.f. = 300/(52 — 2). However, it improves substantiahy if we exclude 
the data of Protopopescu et al., which are strongly biased (as is clearly seen in Fig. 8); the central values, 
however, vary very little. We have included in the errors in (2.19b) half the difference between the two 
possibilities (with or without the data of Protopopescu et al.) 
The fit returns the values 

= (18.0 ± 2.8) X 10-4 x M-^ = (-8.1 ± 3.1) 10"^ x M"^; 
Ff^ = 190 ± 8 MeV . 

The value of 62"^ agrees within 1.3 cr with the more precise result obtained with the Proissart-Gribov repre- 
sentation (see below), that gives 

4°^ = (18.70 ± 0.41) 10"^ X M-^ 6^°^ = (-4.16 ± 0.30) 10"^ x M"^. (2.21) 

This agreement is remarkable, taking into account that the D-wave &2 calculations are less reliable since they 
are the ones with a relatively large contribution of the derivative of the 1=2 t-channel amplitude, which is 
very uncertain (see ref. 6). 

We then take into account the inelasticity by writing 

1, s < 4m^, 

^f(^) = <! l_e-^ , = 0.262 ±0.030; s > Am^- (2-22) 

fc2(M2j' ^' 
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k2 = i/s/4 — m|f. We have fixed the coefficient e fitting the inelasticities of Protopopescu et al.,!^"' and the 
experimental inelasticity of the f2', the overall /d.o.f. of this fit is ~ 1.8. 

In principle, by including b'^^ in our fit the nominal uncertainties in the "fit to data" parameters 
could be decreased. However, we have not done so, since, as explained before and in ref. 6, the D-wave 62 's 
arc less reliable than the 02 's (used as input with very large errors) and, therefore, the apparent improvement 
would be made at the cost of reliability. 

2.4. The F, G waves 



The contribution of the F wave to our sum rules and dispersion relations is very small, but it is interesting 
to evaluate it (we have included its contribution to all relations) to check, precisely, the convergence of the 
partial wave scries. The contribution of the GO, G2 waves is completely negligible, for the quantities we 
calculate in our paper, but we describe fits to these for completeness. 

2.4.1. The F wave 

The experimental situation for the F wave is somewhat confused. According to Protopopescu et al.,f^°l 
it starts negative (but compatible with zero at the 2 a level) and becomes positive around s^/^ = 1 GeV. 
Hyams et al.'^^"' report a positive 63(3) when it differs from zero (above s^^'^ = 1 GeV). In both cases no 
inelasticity is detected up to s^/-^ ~ 1.4 GeV. 

The corresponding scattering length may be calculated with the help of the Proissart-Gribov repre- 
sentation and one finds 

as = (6.3 ± 0.4) X lO""'^ M^^. 

It could in principle be possible that 63(3) changes sign twice, once near threshold and once near s^/^ = 
1 GeV. However, we disregard this possibility and write, simply, 

1/2 

with Sg = 1.45 GeV, and impose the value of 03. 

It is to be understood that this parametrization provides only an empirical representation of the 
available data, and that it may not be reliable except at very low energies, where it is dominated by the 
scattering length, and for 3^^"^ > 1 GeV. We fit data of Protopopescu et al.,[^°l for energies above 1 GeV, 
and data of Hyams et al.I^^^' We have estimated the errors of this last set (not given in the paper) as the 
distance from the average value to the extreme values in the different solutions given. We find 

■^j ~ Yj^, Bo = (1.09 ± 0.03) X 10^ Bi = (1.41 ± 0.04) x 10^ (2.23b) 

The errors have been increased by including as an error the variation that affects the central values when 
using only one of the two sets of data. We do not include separately the effects of the p3(1690), since its tail 
is incorporated in the fitted data. 

2.4.2. The G waves. 

The experimental information on the G waves is very scarce. For the wave G2, we have two nonzero values 
for (54^' from Cohen ct al..['*l and four significant ones from Losty et al.;!^^! they are somewhat incompatible. 
We then fit the data separately, with a one-parameter formula; we write 

If we fit the data of Losty et al. we find B = (—0.56 ± 0.09) x 10®, while from Cohen et al. we get B = 
(-10.2±3.0)xl0®. Fitting both sets together we find B = -9.1x10*^, and a very poor /d.o.f . = 32/(6-1). 
Enlarging the resulting error to cover 6 ct, we obtain our best result, 

cot Sf\s) = ^'^^^^ B, B = (-9.1±3.3) X 10^ (2.24) 
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This formula can only be considered as valid only for a limited range, 0.8 < s^/^ < 1.5 GeV. In 
fact, from the Froissart-Gribov representation it follows that the G2 scattering length is positive. One has 
^ = (4.5 ± 0.2) X 10"^ M'^, while (2.24) would give a negative value. 

For the GO wave, the situation is similar. However, we know of the existence of a very inelastic 
resonance with mass around 2 GeV. An effective value for the imaginary part of the corresponding partial 
wave may be found in Appendix A. 

3. Forward dispersion relations 

We expect that the scattering amplitudes that follow from the phase shifts (and inelasticities) at low energy 
(s^/^ < 1.42 GeV), and the Regge expressions at high energy, will satisfy dispersion relations since they fit 
well the experimental data and are therefore good approximations to the physical scattering amplitudes. 
In the present Section we will check that this is the case, at low energies (,s^/^ < 0.95 GeV), for three 
independent scattering amplitudes (that form a complete set) , which we will conveniently take the following 
t-symmetric or antisymmetric combinations: tt^tt" — > 7r°7r°, 7r°7r+ — > 7r°7r+, and the amplitude 7f = 1, 
corresponding to isospin unity in the t channel. The reason for choosing these amplitudes is that the 
amplitudes for tt^tt" and 7r°7r+ depend only on two isospin states, and have positivity properties: their 
imaginary parts are sums of positive terms. Because of this, the errors are much reduced for them. This is 
easily verified if we compare with the errors in the dispersion relations for tt^tt" and Tr^Tr^ with those for the 
amplitude with It = 1, which depends on three isospin states and has no positivity properties (see below, in 
Fig. 15). 

Here we will not cover the full energy range or try to improve the parameters by requiring fulfillment 
of the dispersion relations, something that we leave for Sect. 4. We will start discussing the global fit in 
Subsect. 2.2.2; the results using the individual fits for SO, as in Subsect. 2.2.3, will be discussed later, in 
Sect. 3.3. 

In our analysis one should take into account that, for the amplitudes that contain the SO wave, the 
uncertainties for it above 0.95 GeV induce large errors in the dispersive integrals, and the agreement between 
dispersive integrals and real parts of the scattering amplitudes evaluated directly becomes affected. This is 
particularly true for the 7r°7r° amplitude, dominated by the SO contribution, where the mismatch becomes 
important above ~ 0.7 GeV. For the 7r°7r+ amplitude, however, since it is not affected by the SO problem, 
the fulfillment of the dispersion relation is good, within reasonably small errors, up to the very region where 
Regge behaviour takes over, s^/^ ~ 1.42 GeV. 

This is a good place to comment on the importance of the contribution of the Regge region (i.e., 
from energy above 1.42 GeV) to the various dispersive integrals. Of course, this depends on each dispersion 
relation. As an indication, we mention that for the unsubtracted dispersion relation (3.7) at threshold, the 
contribution of the Regge region (s > 1.42 GeV) is 9%. For the subtracted relation (3.1a) it is of 10% at 
sV2 = 0.5 GeV and 20% at s^/^ = 0.8 GeV. 

A final general comment is that here, as well as for the sum rules that we will discuss in Sect. 5 and 
the Froissart-Gribov calculations (Sect. 6), we only include the contributions of waves up to and including the 
F wave. We have checked, in a few typical cases, that the contributions of the GO, G2 waves are completely 
negligible. 

3.1. Even amplitude dispersion relations (with the global fit) 
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Figure 9. The combination ReFoo(s) — Foo(4M^) (continuous line) and the dispersive integral (bro- 
ken line). 



3.1.1. TT^TT*' scattering 



We consider here the forward dispersion relation for tt^tt" scattering, subtracted at threshold, 4Af^. We 
write, with Foq{s) the forward tt'^tt*' amplitude, 



ReFoo(.)-Foo(4Aa = ^(i^i^P,^ ^ (2.' - 4M^) Inr ^oo(.') 



s'(s' - s) is' - 4M2) (s' + s - 4Af2) ■ 



(3.1a) 



In particular, for s = 2M^, which will be important when we later discuss the Adler zeros (Sect. 4), we have 



8Mt 



ImFoo(s) 



(3.1b) 



TT y4M^ s(s-2Af2)(s-4M2)- 

We first check the sum rule (3.1b). We take for Foo(4A'/2)^ Foo(2Af2) the values that follow from our 
fits to experimental data of Sect. 2, which provide a representation of partial waves valid below threshold 
(provided s > 0), and evaluate the dispersive integral with the parametrizations we obtained also in Sect. 2. 
We find fulfillment to less than 1 a: 

L»oo (43 ± 3) X IQ-^ (3.2a) 

and 

Foo(4A//2, 0) - Foo(2Af2, Q) = (33 ± 22) x IQ-^. (3.2b) 

For the difference (which should vanish if the dispersion relation was exactly fulfilled, and which takes into 
account correlations) we get 

Foo(4Af2, 0) ~ Foo(2Af2, Q) - Dao = (-10 ± 23) x IQ-^. (3.2c) 

We can also verify the dispersion relation (3.1a). The result is shown in Fig. 9, where a certain 
mismatch is observed in some regions. As we will see below, the matching is better for the other dispersion 
relations because of the smaller weight of the SO wave there. 
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FlGURE 10. The combination ReFo+(s) — _fo+(4Af^) (continuous line) and the dispersive integral 
(broken line). 



3.1.2. vr'^TT^ scattering 

We have, with i^o+(s) the forward tt'^tt^ amplitude, 

Lf^ 5'(.s' - ,)(s' - 4Af2)(s' + s - 4M2) 
In particular, at the point s = SAf^, this becomes 



s(s-2M2)(s-4A.f2)- 

The calculation is now more precise because I3o+ is dominated by the P wave, very well known. We 
find, for the dispersive evaluation, 

130+ = (12± 1) X 10"^. (3.4a) 

On the other hand, using directly the explicit parametrizations for the partial wave amplitudes we found in 
Sect. 2 one has 

Fo+ {AMI , 0) - Fo+ (2 Af 2 , 0) = (6 ± 16) X 10-\ (3.4b) 

For the difference, 

Fo+(4Af2, 0) - Fo+{2Ml,0) - Do+ = (-6 ± 17) x 10"^ (3.4c) 

i.e., perfect agreement. 

This is a good place to remark that the agreement of the values of Fq+{2M'^,{)) and Fqq{2M'^^Q) 
obtained with our parametrizations, and those found evaluating dispersion relations (Eqs. (3.2c), (3.4c)) 
provides a nontrivial test of the validity of our parametrizations even in regions below threshold, well beyond 
the region were we fitted data. 

The fulfillment of the dispersion relation (3.3a) is shown in Fig. 10 for s^/^ below 1.4 GeV. The 
agreement is now good in the whole range; the average x^/j.o.f. for s^/^ < 0.925 GeV is of 1.7. The fact 
that the fulfillment of the dispersion relation reaches the energy where the Regge formulas start being valid, 
s^/2 ^1.4 GeV, is yet another test of the consistency of the Regge analysis with the low energy data. 

To test the dependence of our results on the point at which we effect the junction between the 
phase shift analyses and the Regge formulation, we have repeated the calculation of the forward 7r*'7r+ 
dispersion relation performing this junction at 1.32 GeV, instead of doing so at 1.42 GeV. The fulfillment 
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Figure 11. The amplitude Re_F^* ^(^,0) (continuous line) and the dispersive integral (broken line). 



of the dispersion relation improves slightly (below the percent level) at low energy, s^/^ < 0.75 GeV; while 
it deteriorates a bit more for s^/^ > 0.75 GeV. The net results are practically unchanged, with the choice 
s^^"^ = 1.42 GeV for the junction slightly favoured. We will use this number (1.42 GeV) henceforth. 

3.2. The odd amplitude F^^*—^^; dispersion relation and Olsson sum rule (global fit) 

We consider first a forward dispersion relation for the amplitude with isospin 1 in the t channel, 

evaluated at threshold. This is known at times as the (first) Olsson sum rule. Expressing (4M^, 0) 

in terms of the scattering lengths, this reads 

o (0) r (2) n n r A Im-F(^'=')(g,0) 

2a^ ^ - 5a^ ^ = Doi., Dqi. ^ i^U ds . (3.5) 



4A/2 



s{s - 4 A/2) 



In terms of isospin in the s channel, 

F(^*=i)(s,t) - i/^(^==o)(s,i) + iF(^»=i)(s,t) - lF('^-^Hs,t). 

Substituting in the right hand side above the scattering amplitudes we have just determined up to 
1.42 GeV, and the Regge expression for rho exchange of Appendix B at higher energies, we find, 

Doi. = 0.647 ±0.021. (3.6a) 
(Here, and in all the numbers for scattering lengths and effective ranges, we will take the pion mass as 
unity). This is to be compared with what we find from the values of the a^^ we found in the fits of Sect. 2, 

2^°^ - 54^^ = 0.719 ± 0.072. (3.6b) 

For the difference, 

24°^ - 54^^ - Doi, = 0.073 ± 0.077, (3.6c) 

thus vanishing within errors. 

One can also evaluate the corresponding dispersion relation. 
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evaluating ReF(^'~^)(s, 0) at all values of s, either directly using the fits of Sect. 2, or from the dispersive 
integral. The agreement, as shown in Fig. 11, is very good below 1 GeV, and reasonably good above this. 

3.3. Dispersion relations using the individual fits to SO wave data 

In this Subsection we present the results of checking the forward dispersion relations using the individual 
fits for the SO wave that we performed in Subsect. 2.2.3. The methods are identical to those used for the 
solution with the global fit for this wave in the previous Subsections, so we will skip details and give only 
the results, summarized in Table 2. Here we give the average^ x^/d.o.f. for the dispersion relations for the 
amplitudes that contain the SO wave: It — ^ and tt'^tt'^. The values of the parameters Si, /io are, of course, 
as in Table 1, but we repeat them here for ease of reference. 

We have separated in Table 2 the fits which produce a total x^/d.o.f. of less than 6, from those 
that give a number larger than or equal to 6 when running the corresponding amplitudes through dispersion 
relations. We may consider that the second set is disfavoured by this test. Also, we may repeat some of the 
comments made in Subsect. 2.2.3 with regard to Solution E of Grayer et al., and the evaluation using K 
decay data alone: their errors are very large, due of course to the small number of points they fit, so that 
their fulfillment of dispersion relations is less meaningful than what looks at first sight. 

4. Improving the parameters with the help of dispersion relations 

In this Section we will show how one can improve the results for the fits to the individual waves that we 
found in Sect. 2: the fact that the dispersion relations are fulfilled with reasonable accuracy at low energy, 
and that at the ^ 3 a level they still hold to higher energies, suggests that we may improve the values of the 
parameters we have found with our fits to data requiring also better fulfillment of such dispersion relations. 
This will provide us with a parametrization of the various waves with central values more compatible with 
analyticity and s — u crossing. This method is an alternative to that of the Roy equations to which it is, 
in principle, inferior in that we do not include s — t crossing (although we check it a posteriori); but it is 
clearly superior in that we do not need as input the values of the scattering amplitude for |t| up to 30M^, 
where the various Regge fits existing in the literature disagree strongly one with another (see Appendix B) 
and also in that, with dispersion relations, we can test all energies, whereas the Roy equations are only valid 
for s^/^ < \/60 Af^ 1.1 GeV -and, in practice, only applied up to 0.8 GeV. 

4.1. Improved parameters for the global fit of Subsect. 2.2.2 

We will consider the displacement of the central values of the parameters, requiring fulfillment, within errors, 
of all three dispersion relations for V2M^ < s^/^ < 0.925 GeV (note that we even go below threshold), 
starting with the global solution in Eq. (2.14). We do not fit higher energies because we feel that the errors 
in the input for some waves is too poorly known to give a reliable test there. Specifically, the P wave in 
the region 1.15 GeV < s^^^ < 1.5 GeV is not at all determined by experiment; depending on the fit, a 
resonance appears, or does not appear, in that region: its mass varies between 1.25 and 1.6 GeV. Something 
similar occurs for the SO wave above 0.95 GeV. Thus, it may well be that the parametrizations we use (which, 
for example, assume no P wave resonance below 1.5 GeV) are biased. This could likely be the explanation 
of the slight mismatch of dispersion relations above 1 GeV, when the corresponding contributions begin to 
become important; especially for the rcai parts of the scattering amplitudes. 

Because of this uncertainty with the P wave above 1.15 GeV, and the SO wave above 0.95 GeV, 
which goes beyond their nominal errors, and because they are of small importance at low energy, we have 
not varied the parameters that describe these waves here. We have then minimized the sum of x^'s obtained 
from the variation of the parameters of the waves, within their errors, as obtained from data, plus the average 
of the dispersion relations (that we call "X^ /d.o.f."). This average is obtained evaluating each dispersion 

^ That is to say, the sum of the of each point, spaced at intervals of 25 MeV, divided by the number of points 
minus the number of free parameters. 
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So 


Si 


/io (MeV) 


do.f.(^' = ^) 


2 

A. ^_()_o^ 
d.o.f.^"" ) 


PY, Eq. (2.14) 


21.04 


6.62 


782 ± 24 


0.3 


3.5 


K decay only 


18.5 ± 1.7 


= 


766 ± 95 


0.2 


1.8 


K decay data 
+ Grayer, B 


22.7 ± 1.6 


12.3 ±3.7 


858 ± 15 


1.0 


2.7 


K (lofav (lata 
T (TraytT. (,' 


i().8 ± 0.85 


-().:', i ± 2..:U 


787 ± 9 


0.1 


1.0 


K decay data 
+ Grayer, E 


21.5 ±3.6 


12.5 ± 7.6 


1084 ± 110 


2.1 


0.5 


K decay data 
+ Kaminski 


27.5 ±3.0 


21.5 ±7.4 


789 ± 18 


0.3 


5.0 


K decay data 
+ Grayer, A 


28.1 ± 1.1 


26.4 ± 2.8 


866 ±6 


2.0 


7.9 


K decay data 
+ EM, s— channel 


29.8 ± 1.3 


25.1 ±3.3 


811 ±7 


1.0 


9.1 


K decay data 
+ EM, t-channel 


29.3 ± 1.4 


26.9 ± 3.4 


829 ±6 


1.2 


10.1 


K decay data 
+ Protopopescu, VI 


27.0 ± 1.7 


22.0 ±4.1 


855 ± 10 


1.2 


5.8 


K decay data 
+ Protopopescu, XII 


25.5 ± 1.7 


18.5 ±4.1 


866 ± 14 


1.2 


6.3 


K decay data 
+ Protopopescu, VIII 


27.1 ±2.3 


23.8 ± 5.0 


913 ± 18 


1.8 


4.2 



Errors as in Eq. (2.14b). 

PY, Eq. (2.14): our global fit of Subsect. 2.2.2. The next rows show the fits to K decay'^^' alone or combined with tttt 
scattering data. Grayer A, B, C, E: the solutions in the paper of Grayer et al.'^^"' EM: the solutions of Estabrooks and 
Martin. '^^"^ Kaminski refers to the papers of Kaminski et al.'^^"^' Protopopescu VI, XII and VIII: the corresponding 
solutions in ref. 10. 

Table 2 



relation at intervals of 25 MeV in s^/^, from threshold up to s^/^ = 0.925 MeV, dividing this by the total 
number of points. For the dispersion relations for vr^Tr"*" and tt'^tt^ scattering, we also include in the fit the 
relations (3.1b) and (3.3b), which are important in fixing the location of the Adler zeros for the SO, S2 waves. 

According to this, we allow variation of the parameters of the SO wave up to KK threshold (including 
the location of the Adler zero, zo); the parameters of the P wave up to .s^/^ = 1.0 GcV: and the parameters 
of S2, DO, D2 and F waves for ah s^/^ < 1.42 GeV. For S2 we also leave free Z2. We find that the total 
variation of the parameters has an average of 0.38, showing the remarkable stability of our fits. The only 
parameters that have varied by ^ 1 cr or a bit more arc some of the parameters for the SO and D2 waves. 
For both, this hardly affects the low energy shape, but alters them a little at medium and higher energies 
(for D2, see Fig. 7). Given the low quality of experimental data in the two cases, this feature should not be 
surprising. 

As stated above, in the present Subsection we take as starting point the SO wave we obtained with 
our global fit in Subsect. 2.2.2. The new central values of the parameters, and the scattering length and 
effective range parameters (both in units of M^) are listed below. 
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SO; s^/^ < 2mK : -Bo = 17.4 ± 0.5; Bi = 4.3 ± 1.4; 

Ho = 790 ± 21 MeV; zq = 195 MeV [Fixed]; 
4°^ = 0.230 ± 0.015; 6^°^ = 0.312 ± 0.014. 
S2; < 1.0 : Bq = -80.8 ± 1.7; Bi = -77 ± 5; Z2 = 147 MeV [Fixed]; 
o^^^ = -0.0480 ± 0.0046; 6^^^ = -0.090 ± 0.006. 
S2; 1.0 < < 1.42 : Bo = -125±6; Bi=-119±14; e = 0.17 ±0.12. 

P;s^/2<1.05: = 1.064 ± 0.11; Si = 0.170 ± 0.040; Mp = 773.6 ± 0.9 MeV; 
ai = (38.7 ± 1.0) X 10-2; h = (4.55 ± 0.21) x 10"^ 
DO; s^/^ < . Bo = 23.5 ±0.7; Bi= 24.8 ±1.0; e = 0.262 ± 0.030; 

4°) = (18.4 ± 3.0) X 10"^ 6^°^ = (-8.6 ± 3.4) x 10"''. 
D2; s^/^ < 1.42 : Bq = (2.9 ± 0.2) x 10^; Si = (7.3 ± 0.8) x 10^; B2 = (25.4 ± 3.6) x 10^; 
A = 212 ± 19; 

= (2.4 ± 0.7) X 10"^ b^2 ^ = (-2.5 ± 0.6) x 10~^. 
F; < . Bo = (1.09 ± 0.03) x 10^ Bi = (1.41 ± 0.04) x 10^ as = (7.0 ± 0.8) x 10"^ 
[ 5^/2 < 0.925 GeV; total average xVd.o.f.=0.80 ]. 

(4.1a) 

We note that the central values and errors for all the waves, with the exception of the SO and D2 
waves (and a little the S2 wave), are almost unchanged (in some cases, they are unchanged within our 
precision). The results for the SO, S2 and D2 waves arc shown in Figs. 13, 14 and 7. 

This brings us to the matter of the errors. In general, one cannot improve much the errors we found 
in Sect. 2 by imposing the dispersion relations since, to begin with, they are reasonably well satisfied by our 
original paramctrizations; and, indeed, the errors obtained fitting also the dispersion relations arc almost 
identical to the ones we found in Sect. 2. Any improvement would thus be purely nominal and would be 
marred by the strong correlations among the parameters of the various waves that would be introduced. 

The exceptions arc, as stated, the SO, S2 and D2 waves. For the first two when including the 
fulfillment of the dispersion relations in the fits, we have first left the Adler zeros, located at ^Zq and 2z2, 
free. We find 

zo = 195 ± 21 MeV, Z2 = 147 ± 7 MeV . (4.1b) 

Unfortunately, the parameters are now strongly correlated so the small gain obtained would be offset by the 
complications of dealing with many correlated errors.^ Because of this we have fixed the Adler zeros at their 
central values as given in (4.1b) when evaluating the errors for the other parameters. Then the errors are 
almost uncorrelated. 

For the D2 wave we accept the new errors because its parameters vary by more than 1 a from those 
of Sect. 2. Given the poor quality of experimental data, obvious from a look at Fig. 7, we feel justified in 
trusting more the central values and errors buttressed by fulfillment of dispersion relations. 

We consider (4.1) to provide the best central values (and some improved errors) for the parameters 
shown there. The various dispersion relations are fulfilled, up to = 0.925 GeV, with an average of 0.66 
(for tt'^tt'^), of 1.62 for the 7r°7T+ dispersion relation, and of 0.40 for the /« = 1 case. The consistency of our tttt 
scattering amplitudes that this shows is remarkable, and may be seen depicted graphically for the dispersion 
relations in Fig. 12, where we show the fulfillment of the dispersion relation even up to s^/^ = 1.42 GeV. 
Mismatch occurs to more than one unity of x^/d.o.f . due to the artificial joining of our low and higher energy 
fits to data, and, above all, to the uncertainties of the P and, especially, the SO wave between 1 and 1.42 
GeV. 

* There is another reason for not attaching errors to zo; the Adler zero is located at ^Zq ~ 0.01 GeV^, so near the 
left hand cut that our conformal expansion cannot be considered to be convergent there. 
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Figure 12. Fulfillment of dispersion relations, with the central parameters in (4.1a). The error bands 
are also shown. 
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In fact, this mismatch is small; if we take the improved values of the parameters as given in (4.1a), 
and rc-calculatc the various dispersion relations up to s^/^ = 1.42 GcV, wc find that for tt^tt^ scattering the 
dispersion relation is fulfilled with an average of 1.85, while for tt'^tt^ scattering we find an average of 
1.57 and the /f = 1 dispersion relation is fulfilled to an average of 1.16. However, the last two numbers 
become smaller than unity if only wc increase the error of the slope of the P wave between 1 and 1.4 GeV 
by a factor of 2, i.e., if we take Ai = 1.1 ± 0.4 in Eq. (2.7). 

Likewise, if we replaced the SO inelasticity in (2.15b) by that in (2.15b'), the average for tt^tt^ 
would improve to 1.35, while the average for the It = I amplitude would become slightly worse, 1.47. 
Doubtlessly, the incompatibility of phase and inelasticity for the SO wave above 1 GeV precludes a perfect 
fit. We plan, in a coming paper, to use the results here to find consistent values for the phase shift and 
inelasticity for the SO and P waves between 1 GcV and 1.42 GeV, by requiring consistency below the la 
level of the dispersion relations in that energy range. 

The agreement of dispersion relations up to 1.42 GeV is the more remarkable in that the dispersion 
relations above 0.925 GcV have not been used to improve any wave. 

For the sum rules (3.2) and (3.4) we now find 

24°^ - 54^^ - Doi. = (25 ± 32) X 10"^ 

Foo(4M2,0) -Foo(2M2,0)--Doo = (-15±9)x10-3, (4.2) 
i^o+(4M^,0) - Fo+{2MlO) - Do+ = (3± 7) x lO'^. 

4.2. Improved parameters for the individual fits of Subsect. 2.2.3 

We now refine the parameters of the fits to data, but taking as starting point the numbers obtained from 
the individual fits to the various data sets as described in Subsect. 2.2.3. Apart from this, the procedure is 
identical to that used in Subsect. 4.1; in particular, we also fit the relations (3.1b) and (3.3b). The results 
of the evaluations are presented in Table 3, where we include the solution (4.1), and also what we find if 
requiring the fit to only K decay data for the SO, as given in Eq. (2.12). 



Improved fits: 


Bo 


Si 


(McV) 


zo (McV) 


x'ih = 1) 
d.o.f. 


2/ 0\ 

X (vr TT ) 
d.o.f. 


d.o.f. 


(3.1b) 


FY, Eq. (4.1) 


17.4 ±0.5 


4.3 ± 1.4 


790 ± 21 


198 ± 21 


0.40 


0.66 


1.62 


1.6 a 


K decay only 


16.4 ±0.9 


= 


809 ± 53 


182 ± 34 


0.30 


0.29 


1.77 


1.5 a 


K decay data 
+ Grayer, C 


i(i.2 ± 0.7 


0.5 = 1.8 


788 ± 9 


181 ± :!<) 


().:!7 


().:i2 


1.71 


1. ") a 


K decay data 
+ Grayer, B 


20.7 ± 1.0 


11.6 ±2.6 


861 ± 14 


233 ± 30 


0.37 


0.83 


1.6 


4.0(7 


K decay data 

+ Grayer, E 


20.2 ±2.2 


8.4 ±5.2 


982 ± 95 


272 ± 50 


0.60 


0.09 


1.4 


6.0(7 


K decay data 
+ Kaminski 


20.8 ± 1.4 


13.6 ±3.7 


798 ± 17 


245 ± 39 


0.43 


1.08 


1.36 


4.5 (T 



FY, Eq. (4.1): our global fit, improved with forward dispersion relations, as explained in Subsect. 4.1. Grayer A, B, 
C, E means that we take, as experimental low energy data for the SO wave, the solutions in Grayer et al.,'^^"' fitted 
as shown in Table 1. Kaminski means wc have used the data of Kaminski ct al.'^^''' In these fits we have improved 
the parameters requiring fulfillment of dispersion relations. Although errors are given for the Adler zero, we have 
fixed it at its central value when evaluating other errors. We have included the fulfillment of the sum rule (3.1b) in 
the last column, but not (3.3b) and (3.5), which are verified within 1 cr by all solutions. 

Table 3 
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FIGURE 13. The SO phase shift corresponding to Eq. (4.1) (PY, thick continuous line and error 
band), the unimproved solution of Eq. (2.14) (thin continuous line), and the improved solutions "K 
decay only" and "Grayer C" of Tabic 3 (difficult to see as they fall almost on top of PY). The solution 
CGL'^' (lowest discontinuous line) is also shown. 

The following comments are in order. First of all, we have the remarkable convergence of the first 
three solutions in Table 3: and even the last three solutions approach our evaluation, PY, Eq. (4.1). This 
convergence is not limited to the SO wave: the parameters for all waves other than the SO agree, within 
< 1(7, for all solutions in Table 3, with the numbers given for the PY solution in Eq. (4.1). For this reason 
we do not give in detail, for each individual fit, the improved solutions for waves other than the SO wave. 
Of course, the values of the x^'s given in Table 3 have been evaluated with the corresponding improved 
solutions for all waves. 

Secondly, the fulfillment of dispersion relations is similar for all solutions in Table 3. However, the 
solutions based on Grayer B, Grayer E, and that based on the data of Kaminski et al., fail to satisfy the 
sum rule (3.1b) by a large amount even though they have large errors; sec the last column in Table 3. We 
should, therefore, consider the three solutions PY, that based on K decay data only, and the solution based 
on Grayer C, to be clearly favoured by this consistency test. This is satisfactory in that the solution C was 
obtained by Grayer et al.I"^"'^"! from solution B by including absorption corrections. 

Thirdly, we may consider the solution based on K decay data only, and with a single parameter Bq, 
to be in the nature of a first approximation, and the introduction of the parameter Bi as producing the more 
accurate solutions denoted by Grayer C and PY. In this respect, the coincidence of the common parameters 
of the first three solutions in Table 3 within < 2 cr is satisfactory. And, indeed, this coincidence is even 
more pronounced. In Fig. 13 we show our starting, global fit from Eq. (2.14) together with the improved 
solution PY, Eq. (4.1) and the improved solution corresponding to Grayer C: they arc all three contained 
inside the error band of PY, Eq. (4.1), with which the solution Grayer C practically overlaps. We also show, 
in Fig. 14, the S2 wave, before and after improvement. 

Finally, and in spite of the virtual coincidence between the two improved solutions Grayer C and 
PY, Eq. (4.1), we consider the last to be preferred: it incorporates data from various experiments and is thus 
less likely to be biased by systematic errors. Because of all this, we feel confident in considering our solution 
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FlGURE 14. Phase shift for the S2 wave, Eq. (4.1) (PY, thick continuous line, and error band), our 
unimproved fit, Eq. (2.8) (dotted line) and the solution CGL,'^' and error band (dashed line). We also 
show the data from Durusoy et al.'^^' and solution B in Hoogland et al.,'^^' not included in the fit. 

in Eq. (4.1) to be fully validated, and the method we have used to be well tested. We will then henceforth 
accept and work with the solution given in Eq. (4.1) for the improved SO wave below 0.95 GeV. 

5. Consistency tests of s — t — u crossing: two sum rules 



In this Section we discuss two sum rules that follow from crossing symmetry; in next Section we will consider 
Froissart-Gribov sum rules, which can also be viewed as checks of s — t — u crossing. The sum rules we 
discuss now relate high (s^/^ > 1.42 GeV) and low energy, with the low energy given by the P, D, F waves 
in the region s^^^ < 1.42 GeV and the high energy is dominated by, respectively, the rho and Pomeron 
Rcgge trajectories. The interest in checking them is that they were given by the authors in refs. 1, 2 
(following Penningtonl-'^^l ) as the reason for the incorrect Regge parameters. Contrarily to the assertions in 
these references, however, we will here check, once again, that there is perfect consistency provided one uses 
standard Regge behaviour for energies above 1.42 GeV, and accurate representations of the experimental 
partial waves below 1.42 GeV, such as the ones found here in Sects. 2, 4 and collected in Appendix A. We 
should perhaps remark here that the contributions of the SO, S2 waves cancel in both sum rules, hence we 
do not even need to worry about which solution to use for the SO wave. 

The first sum rule is obtained by profiting from the threshold behaviour to write an unsubtracted 
forward dispersion relation for the quantity F'-^=^^\s, 0)/(s — 4Af^). One gets 

6AU _ 1 , lmf(^-=i)(s,0) , l^^(su) Imf-W(s,0) 



which is known at times as the (second) Olsson sum rule; see e.g. the textbook of Martin et al.^ The index 
1 refers to isospin in the s channel and C^^"'' are the s ~ u crossing matrix elements. Canceling a\ with the 



^ Martin, B. R., Morgan, D., and and Shaw, G. Pion-Pion Interactions in Particle Physics, Academic Press, New York 
(1976). 
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Proissart-Gribov expression for this quantity (cf. Sect. 6), and substituting the c[f'\ we find the result 

f-^^In.F«--)(MM;)-ImFl'-)(..0)_ f d/ff" '^^l Imf «---"(,.0) = 0. (5.2) 

The contributions of the S waves cancel in (5.2), so only the P, D, F and G waves contribute. At high 

energy the integrals are dominated by rho exchange. We take the improved central values of the parameters 
for the different waves as given in Eq. (4.1). In units of the pion mass we get 

(-0.12 ± 1.27) X 10-^ (5.3) 

that is to say, perfect consistency. 

The second sum rule we discuss is that given in Eqs. (B.6), (B.7) of ref. 1. It reads. 



Jami y 



Here F'^^\s,t) = dF^^\s,t)/d cos 9, and the upper indices refer to isospin in the s channel. We get, again 
with the improved parameters and with M,r = 1, 

J= (-0.2±4.2) X 10"^ (5.5) 

These are not the only crossing sum rules that our tttt amplitude verifies; the coincidence of the values 
for the parameters ai, bi, b'^ ^ obtained from direct fits to data in Sect. 2 with those from the Froissart- 
Gribov projection, that involves simultaneously s, u and t crossing, are highly nontrivial ones. We will see 
this in next Section. 

6. Low energy parameters for waves with / > 1 from the Froissart— Gribov 
projection, and a new sum rule for bi 

6.1. The Froissart— Gribov representation for ai, bi 

The quantities oi , bi may be evaluated in terms of the It = 1 amplitude using the Froissart-Gribov repre- 
sentation (for more details on the Froissart-Gribov representation, see refs. 6, 16 and work quoted there). 
For, e.g. ai with I =odd, we have 

o^(/=i)_ V^ni + ^) r ,„ ImF(^*=^)(.,4M^) 



AM^rg + 3/2) . 



r d.'"^'""''''""-'. (6.1) 



We will use the parametrizations in (4.1) for the imaginary parts of the scattering amplitudes in all the 
Froissart Gribov integrals, and the Regge expressions in Appendix B here at high energy. We find the 
following results, in units of M^: 

F.-G. Eq. (4.1) 
lO^xai: 37.7 ± 1.3 38.7 ±1.0. ^ ' 

We here compare the result obtained with (6.1), denoted by F.-G., and the value for ai which we found in 
our improved fits, as given in Eq. (4.1a). Because the two determinations are essentially independent, we 
can compose the errors to get a precise and reliable value for ai : 

ai = (38.4 ± 0.8) x 10"^ M'^. (6.3) 

For the quantity bi we have, with the same conventions as before, 

F.-G. Eq. (4.1) 
10^x61: 4.69 ±0.98 4.55 ±0.21. ^ 

The Froissart-Gribov integral here is dominated by the rho Regge pole, as the low energy contributions 
cancel almost completely. The agreement between these two determinations of 61 is, therefore, a highly 
nontrivial test of the consistency of the high and low energy parts of our pion-pion scattering amplitude. 
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unfortunately not very precise because of the large error in the Froissart-Gribov determination. We will 
discuss this further at the end of Appendix B. 

The F wave scattering length may be similarly evaluated; we find 

as = (6.3 ± 0.4) X 10"^ M"^. (6.5) 



6.2. The Froissart-Gribov projection for even amplitudes: the j 62 parameters 



We first calculate the two combinations of scattering lengths a(0+) = §[02 "^2 ] anda(OO) = |[a2 +2a2 ]• 
They correspond to the s— channel amplitudes 

for which (as mentioned before) errors are minimized. 

The dominant high energy part in the Froissart-Gribov representation is given now by the Pomer- 
anchuk trajectory and its importance is smaller than previously because the integrals converge faster. We 
find, in units of 

a(0+) = (10.61 ± 0.14) X 10"^ M'^- a(00) = (16.17 ± 0.75) x 10"* M'^. (6.6) 
For the effective range parameters, 

6(0+) = (-0.183 ± 0.061) X 10"^ M"^; 6(00) = (-7.96 ± 0.57) x lO"'' M"^. (6.7) 



6.3. A new sum rule for 61 



The Froissart-Gribov representation for the effective range bi depends strongly on the Rcggc parameters 
for rho exchange, and is affected by large errors. Here we will present a sum rule that, contrarily, depends 
almost entirely on the low energy scattering amplitudes and is much more precise. It is obtained in a way 
similar to that used for the first crossing sum rule in Section 5. We now write a dispersion relation for the 
quantity 

d_ / f(^°=i)(s,o) \ 

ds\s- AMI J ' 
which we evaluate at threshold. Taking into account that 

d_ / F(^-=i)(5,0) \ ^ 3M^ 

ds \ 8 -AMI )s=AMl 27r 

we obtain a fastly convergent relation for h\: 



M^bi 



J 4, 



ds 



4M2 
1 



(S - 4M2)3 s3 



ImF(^*=o)(s,0) + | 



+ 



(S-4M2)3 S' 



ImF(^*=^)(s,0) 



(s - 4M2)3 s3 



ImF(^*=2)(s,0) i. 



(6.8) 

Most of the contribution to bi comes from the SO and P waves at low energy, while all other contributions 
(in particular, the Regge contributions) are substantially smaller than 10~^. Adding all pieces we find 

61 = (4.99 ± 0.21) X 10-3 M-^ (6.9) 

a value reasonably compatible with what we found in (4.1), bi = (4.55 ± 0.21) x 10"^ (because of 

correlations, the distance is actually below 1 a). We can combine both and find a precise estimate, 

bi = (4.75 ± 0.16) X 10-3 M-^. (6.10) 
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7. Comparison with the low energy parameters of CGL, DFGS and KLL 

We present in Table 4 a global comparison of the low energy parameters as given here, that we denote 
by PY, as well as recent evaluations, that use the Roy equations, by Colangelo, Gasser and Leutwylerl^l 
(CGL), by Descotes at alJ^l, that we denote by DFGS, and by Kaminski, Lesniak and Loiseaul^l, denoted 
by KLL. Besides scattering lengths and effective range parameters, we give the quantities [oq"^ — Oq^]^ and 
(5o°^(m|-) — 5^\m\), relevant, respectively, for pionic atom decays and CP violating kaon decays. 
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17.5 ±0.3 


18.70 ±0.41 (xlO-*) 


af 






1.70 ±0.13 


2.78 ±0.37 (xlO"") 


a(0+) 








10.61 ±0.14 (xlO^*) 


a(00) 








16.17 ±0.75 (xlO^*) 


bf' 






-3.55 ±0.14 


-4.16 ±0.30 (xlO-*) 


6f 






-3.26 ±0.12 


-3.89 ±0.28 (xlO-*) 


6(0+) 








-0.183 ±0.061 (xlO-*) 


6(00) 








-7.96 ±0.57 (xlO""*) 


as 






5.6 ±0.2 


6.3 ±0.4 (xlQ-^) 



Units of Mtt- The numbers in the CGL column are as given by CGL in Table 2 and elsewhere in their text. 

In PY, the values for the D, F waves parameters arc from the Froissart Gribov representation. The rest are from the 
fits, improved with dispersion relations, except for ai and 61 that have been taken from Eqs. (6.3) and (6.10). 

Table 4 
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Figure 15. Consistency of disper- 
sion relations for the tttt amplitudes of 
ref. 2 (CGL) and for our amplitudes, 
with the parameters in (4.1), denoted by 
PY. We plot the differences Ai, given in 
Eqs. (7.1), between the results of the cal- 
culation of the real parts directly with 
the various parametrizations given, or 
from the dispersive formulas. Perfect 
consistency would occur if the continuous 
curves coincided with the dotted lines. 
The error bands are also shown. The 
progressive deterioration of the CGL re- 
sults as the energy increases is apparent 
here. 



The DFGS solution is compatible, within its errors, both with CGL and PY. As to the KLL solution, 
its values for a'^^\ bg^^ are compatible with those of DFGS (for the first), somewhat less so with the numbers 
of PY and with what CGL find. The central value for ai of KLL is too high, although it is compatible 
within its errors with other determinations. The value of 6i, however, is Str below the next lowest one in 
the Table 4. The reason could be that Kamiiiski, Lesniak and Loiseau use some approximate calculation 
techniques, like effective, separable potentials and Pade approximants. 

The calculation of CGL^^l is different from the others, since CGL impose ch.p.t. to two loops. 
However, CGL use Roy equations, which require phenomcnological input for energies above 0.8 GeV. As 
stated in the Introduction (and discussed in more detail in refs. 6, 7) some of this input is not accurate, 
which may cause biases in the final results. In this respect, wc find that the predictions of CGL for the 
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scattering lengths ag \ ai agree very well with what we find from experiment. This means that our results 
do not challenge the standard chiral counting of ch.p.t. However, the solution given by CGL deviates from 
our results for quantities like 61, ai'\ h^P or the SO and S2 phase shifts above 500 MeV (and in particular 
the value used as input by CGL at the matching point 4°'(800Mey) = 82.3 ± 3.4°). 

That the CGL solution deteriorates as the energy increases is quite transparent if we compare the 
fulfillment of dispersion relations with the parameters of CGL for the SO, S2 and P waves al low energy, or 
with our parameters here. This is depicted in Fig. 15, where we show the mismatch between the real part 
and the dispersive evaluations, that is to say, the differences Ai, 

^1 =ReF(-^'-^) s>0 2LP.P. / ds'— — ^ 1,0^ , 7.1a 

^'^ TT hui (s'-s)(s' + s-4M2)' ^ ' 

and 

Ao+ = ReFo+is)-FoA^M^) P-P-^/^^^d. 4^2) + 4^2) " (^.Ic) 

These quantities would vanish, Ai = 0, if the dispersion relations were exactly satisfied. 

We include in the comparison of Fig. 14 the errors; in the case of CGL, these errors are as follow 
from the parametrizations given by these authors in ref. 2, for s^/^ < 0.8 GeV. At higher energies they 
are taken from experiment via our parametrizations. By comparison, we show the same quantities for our 
best results in the present paper, that is to say, with amplitudes improved by use of dispersion relations, 
Eq. (4.1). In both cases we have taken the Regge parameters from Appendix B here. 

8. Summary of conclusions 



In the previous Sections we have given a representation of the tttt scattering amplitudes obtained fitting 
experimental data below 1.42 GeV, supplemented by standard Regge formulas above this energy. We have 
shown that our representations satisfy reasonably well forward dispersion relations at low energy, as well 
as crossing sum rules. Wc have shown that requiring fulfillment of the dispersion relations up to s^^"^ = 
0.925 GeV leads to a refinement of the central values and errors for the parameters for the various waves, 
giving a set of these quantities such that dispersion relations are now satisfied, within acceptable errors, at 
all energies. 

In particular, for the SO wave, we have analyzed the results found starting from different sets of data. 
If we eliminate from those sets the ones that are less consistent with the dispersion relations, it is seen that the 
remaining determinations converge to a solution, essentially unique, when improved by requiring fulfillment 
of the dispersion relations. We have, therefore, obtained a complete set of tttt scattering amplitudes that are 
consistent, with theoretical requirements as well as with experiment: they are collected in Appendix A. 

After this, we use these scattering amplitudes to evaluate low energy parameters for P, DO, D2, and 
F waves in a reliable manner, clearly improving on previous work.I^'^'^l These parameters may then be used 
to test chiral perturbation theory to one and two loops, or to find quantities relevant for pionium decays or 
CP violating kaon decays: 

- o-o^? = 0.077 ± 0.008 M-2, S^°\m%) - S^^\ml) = 52.9 ± 1.6°. 

We may remark here that our errors are, typically, a factor 2 to 3 times larger than those in CGL,[^1 
at low energy. Given that the scattering amplitudes of this reference show mismatches between high and 
low energy, for many observables by several standard deviations (cf., for example. Fig. 14), and that our 
results arc of higher precision at intermediate energies for some quantities (compare, for example, the error 
bars in Figs. 15), we do not consider our results to be inferior. But one certainly can ask the question 
whether it would be possible to improve on our precision. The answer is, no in the sense that our amplitudes 
agree, within errors, with theoretical requirements and with data. A sizable improvement would require 



- 33 - 



-J. R. PELAEZ AND F. J. YNDURAIN- 



substantially improved experimental data, certainly for the SO, S2 waves at low energy, and for the SO and 

P waves above ^ 1 GcV. One may think that imposing chiral perturbation theory could lead to decreasing 
the errors of the tttt scattering amplitudes. However, the matter is complicated and will be left for a future 
publication. 



-34- 



-THE PION-PION SCATTERING AMPLITUDE- 



Appendix A: Summary of low energy (sV2 < 1.42 GeV) partial waves 

In this Appwidix wc collect the best values for the parainetrizations of the various partial waves. We give the 
values obtained after improving with the help of dispersion relations, Eq. (4.1). The values of the parameters 
and errors obtained from raw fits to data are as given in Sect. 2. 

A.l. The S wave with isospin zero below 950 MeV 

We impose the Adler zero at s = \zq, with zq fixed, and a "resonance" with mass hq, a free parameter. 
After improving the global fit with dispersion relations we have 



cot 5^ '{s) = — — — < Bo + Bi—^ — )■ ; 



(Al) 



Bo=17.4±0.5, Si=4.3±1.4, /io = 790 ± 21 MeV; zq = 195 ± 21 McV [Fixed]. 



This fit we take to be valid for s^/^ < o.95 GeV. Note that we have fixed the value of zq when evaluating 
the errors of the other parameters. 
This gives 

4°^ = (0.230 ± 0.015) M-i; 6^°^ = (0.312 ± 0.014) M'^; S^°\mK) = 44.4° ± 1.5° . (A.2) 
A.2. The J = S wave between 950 MeV and 1420 MeV 
We write 



(U) , I h k'i\ M'"^ - s ,. ^, 

^^^ = l-(^™+^2f ^— . (A.4) 



Then, 



(A.5) 



Co = 1.3 ± 0.5, Mf = 1320 ± 50 MeV; M, = 920 MeV (fixed); 

ei = 6.4 ± 0.5, 62 = -16.8 ± 1.6; M' = 1500 MeV (fixed). 

It should be clear that this fit is not very reliable since there are several incompatible sets of experimental 
data. Thus, the errors arc purely nominal. If using tttt — > KK data to fit the inelasticity we would have 
obtained different values for the ef. 

ei=2.4±0.2, £2 = -5.5 ±0.8. (A.5') 

This fit has not been improved with dispersion relations. For dispersion relations below 0.925 GeV, it is 
irrelevant whether we use (A.5) or (A.5'). 

A.3. Parametrization of the S wave for 1 = 2 below 1 GeV 



We fit only the low energy data, s^/^ < 1.0 GeV, taking the inelastic threshold at sq = 1.05^ GeV 



2 



cot (.) = ^ Ibq + B,^'- ^ ' \ , .V2<i.oGeV. (A.6) 



(A7) 



2k s — 2z2 \ \fs + v'so — s J ' 

The central values and errors are improved with dispersion relations; we get 
= _ 80.8 ± 1.7, Bi=-77±5, Z2 = 147±7MeV [Fixed]; 
af^ =(-0.0480 ±0.0046) x M'^; bf^ = (-0.090 ± 0.006) x M'^; 5^^\mji) = -8.5 ±0.3°. 
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A.4. Parametrization of the S wave for 7 = 2 above 1 GeV 



At high energy, we may fit data from s^l"^ = 0.91 GeV up to s^l"^ = 1.42 GeV. We require junction with the 
low energy fit at 1 GeV, neglect inelasticity below 1.45 GeV, and write 

cot = f^^! J% ji^o + i?r^-^| ; 

° ^ ^ 2fc s - 2^1 \ " ' + J (A.8a) 

3^/2 ^ 45 Qgy. ^ (fixed). 
After improving with dispersion relations we then find 

Bo = -125 ±6, Bi = -119±14; s > (1.0 GeV)^ (A.8b) 
The inelasticity may be described by the empirical fit 

r]P{s) = l-e{l-s/sf/'^, e = 0.17 ±0.12 (5^2 = 1.05 GeV). (A.8c) 
These formulas are expected to hold from 1 GeV up to 1.42 GeV. 

A. 5. The P wave below 1 GeV 

We have 

cot5i(s) = ^(M^ - s) jso + Bi^^^^l ; s^/' = 1.05 GeV. (A.Qa) 
The best result, from (4.1), is 



(A.9b) 



Bo =1.064 ±0.011, Bi = 0.170 ± 0.040, Mp = 773.6 ± 0.9; 
ai = (38.7 ± 1.0) X IQ-^M-^, h = (4.55 ± 0.21) x lO'^M'^ 
Slightly better values for ai, 61, improved with sum rules, may be found in Table 1. 

A.6. The P wave for 1 GeV < s^/^ < 1.42 GeV 

For the P wave between 1 GcV and 1.42 GcV we use an empirical formula, obtained with a linear fit to the 
phase and inelasticity of Protopopescu ct al.l^"l and Hyams ct al.;!^^! 

5i(s)=Ao + Ai(v/^-l), mis) = I - e{^7/§ - 1); 
e = (0.30 ±0.15), Ao = 2.69 ±0.01, Ai = l.l±0.2. 

Here s = 1 GeV^. The value of Aq ensures the agreement of the phase with the value given by (A. 9) at 
s = s = 1 GeV^. This fit is good; however, it should be remarked that there are other sets of experimental 
data for this wave disagreeing with the one used here for s^/^ > 1.15 GeV. Hence (A. 10) may be biased 
beyond its nominal errors. 

This fit has not been improved with dispersion relations. 
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A.7. Parametrization of the D wave for 7 = 



We fit data on 62'^ altogether neglecting inelasticity, which we will then add by hand. The data are scanty, 
and of poor quality. We impose in the fit the scattering length, as obtained from the Proissart-Gribov 
representation, and the experimental width of the f2- Moreover, we improve the central values with dispersion 
relations. 

We write 

1/2 

cot^f (s) = 1^ {Ml - s) + Biw{s)], (A.lla) 

and 

w{s) = V ^o~^ ^ sy^ = 1450 MeV; M/, = 1275.4 MeV [Fixed]. 

We find, 

So =23.5 ±0.7, Bi= 24.8 ±1.0. (A.llb) 
We take into account the inelasticity by writing 

r 1, s< Am\; 

l_e-^, e = 0.262 ±0.030, a > 4m^. (^-H^^) 

Here = -\/s/4 — m^. 

The fit returns the values 

4°^ = (18.4 ± 3.0) X 10"^ X M-^ = (-8.6 ± 3.4) x 10"^ x M"^. 

Better values for b^*^ are given in Table 1. 
A. 8. Parametrization of the D wave for 1 = 2 

For isospin equal 2. there arc no resonances in the D wave. If wc want a parametrization that applies down 
to threshold, we must incorporate the zero of the corresponding phase shift. So we write 

with A a free parameter and 



cotSi'\s) = '^{Bo + BMs)+B2w{sf} ^1 ' _ ' ^'^^ < 1-05 GeV. (A.12a) 



Ks)= ^~^^!°^ , 4/^ = 1450 MeV. 



Moreover, we impose the value for the scattering length that follows from the Froissart-Gribov representation. 
We assume negligible inelasticity and find, after improvement with dispersion relations, 

Bo = (2.9±0.2) X 10^, Bi = (7.3±0.8) X 10^ B2 = (25.4 ± 3.6) x 10^ Zi = 212 ± 19 McV . (A.12b) 

(2) 

The fit returns reasonable values for the scattering length and effective range parameter, 62 • 

4^^ = (2.4 ± 0.7) X 10"^ = (-2.5 ± 0.6) x lO""^ M^-\ (A.13) 

Better values for \ are given in Table 1. 
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A.9. The F wave 

For the F wave below s^/^ = 1.42 GeV we fit the phase shifts of Protopopescu et al.,[^°^ and those of Hyams 
et al.,[^^''l plus the scattering length as given by the Proissart-Gribov representation. After honing the central 
values of the parameters with the help of dispersion relations we find 

Bo = (1.09 ± 0.03) X 10^ Bi = (1.41 ± 0.04) x 10'\ 

This implies 03 = (7.0±0.8) x 10~^ , in agreement (within errors) with the result of the Froissart-Gribov 
calculation (Table 4). 

The contribution of the F wave to all our sum rules is very small (but not always negligible); the 
interest of calculating it lies in that it provides a test (by its very smallness) of the convergence of the partial 
wave expansions. 

A. 10. The G waves 



For the GO wave, wc take its imaginary part to be given by the tail of the /4(2050) resonance, with its 
properties as given in the Particle Data Tables: 

Im f(°)(s) = ( BR • 

[HMI)) '^{s-MlY+Min[k{s)/k{Ml)Y^' (A.15) 

> 1 GcV; BR =0.17 ±0.02, M/^ = 2025 ± 8 MeV, T = 194 ± 13 MeV . 
For the wave G2, we can write, neglecting its eventual inelasticity, 

cot 5f\s) = ^^^y^ B, B = (-9.1±3.3) X lO'^; s^/^ > 1 GeV. (A.16) 

It should be noted that the expressions for the GO, G2 waves, are little more than order of magnitude 
estimates. Moreover, at low energies they certainly fail; below 1 GeV, expressions in terms of the scattering 
length approximation, with 

4°^ = (8.0 ± 0.4) X 10-^ M-^ af ^ = (4.5 ± 0.2) x 10"^ M'^, 

are more appropriate. If, in a given calculation, the contribution of either of the two G waves is important, 
it means that the calculation will have a large error. We have checked in a representative number of our 
calculations that the contributions of the GO, G2 waves are totally negligible. 

Appendix B: The Regge amplitude (s^/^ > 1.42 GeV) 



In the calculations of the various sum rules and dispersion relations we require the imaginary part of the 

scattering amplitudes at high energy. For this we use Regge formulas, applicable when s ^ A^, with 
A ~ 0.3 GeV the QCD parameter, and s :$> \t\. Specifically, we use them above s^/^ = 1.42 GeV and 
expect them to be valid at least for \t\ ^ 4M^ ~ 0.08 GeV"^. These formulas are obtained from A^'A'', wN 
scattering, using factorization, and from direct fits to tttt data. Wc will not include refinements that take into 
account the logarithmic increase of total cross sections at ultra-high energies (that may be found in ref. 5), 
superfluous for applications to tttt scattering. We give here the relevant formulas, for ease of reference, for 
Pomeron and P' Regge poles, and will only discuss in some detail (and even improve slightly) the parameters 
of the rho residue. 
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We consider three Regge poles: Pomeron (P), the P' (associated with the /2(1270) resonance) and 
the rho. For the first two, 

t fixed 

P{s,t) =/3pap(t) i±|^c''*(s/s)"-(*), (B.l) 

api {0)[1 + Q:p'(0)J 

For exchange of unit isospin, we write 

ImFii*=i)(s,i) ^c^^ p{s,t); pis,t) = (B.2a) 

t fixed 

In the present paper we have taken a formula shghtly different from that in ref. 5, viz., 

'^'^^^ = ^'^^^ tTm§ + tSo + ' ^p = ^ + dp + TT^y ^^-^^^ 

In ref. 5, where we took an expression from Rarita et al.^^^'l solution la, assuming the t dependence to be 
equal for tttt and irN scattering (see below, Eq. (B.7)). Our expression now, (B.2b), is simpler. As we will 
show below, we can fix its parameters (for small values of t) without assuming the equality of tttt and TrAi". 
For exchange of isospin two, 

ImPi^'=2)(s,t) R2{s,t) = /32e''*(s/s)"^(*)+"''(°^-^ (B.3) 

t fixed 

This last amplitude corresponds to double rho exchange; the expression we use for it at t ^ is rather 
arbitrary, since little is known about it. Fortunately, it has almost no influence in the calculations. 
Following the analysis of ref. 5, the trajectories are taken as follows: 

ap{t) ~ ap(0) + a'pt, ap{t) ~ ap{0) + a' t+ ^a" t\ (B.4) 

and one has 



-4 



ap(0) =0.52 ±0.02, a^ = 0.90GeV"^ a;;,' = -0.3 GeV" , 
ap(0)=l, ap = (0.20 ±0.10) GcV"^. ''^'^^ 
Moreover, in the present paper, wc have taken the following values for the parameters of the residues: 

/3p =2.54 ±0.03, /?p' = 1.05 ± 0.02, /32 = 0.2 ± 0.2, 6 = (2.4 ± 0.2) GeV"^; 

/3p =/3p(0) = 1.02 ±0.11, dp = (2.4 ± 0.9) GeV"^ ^^'^^ 

The parameters for P, P' are like in ref. 5. Those for the rho trajectory are as improved below. 

The Pomeron and P' parameters arc well fixed from factorization and (besides tttt cross sections) by 
TriV, NN data, down to kinetic energies c:^ 1 GeV; but some extra words will be added on the rho residue, 
Pp{t). This quantity was obtained in ref. 5 in a mixed manner: /3p = /3p(0) was found fitting high energy tttt 
data, improving the error with the help of the sum rule (5.2). The t dependence of (3p{t) was taken as in 
solution la of Rarita et al.,[^''l assuming it to be identical for tttt and ttA'' scattering. Thus, we wrote 

m = /3p(0) [(1.5 + l)e''* - 1.5] 112m t^o + ' = 2-5 + YT^Y ^^'^^ 

This, in fact, is similar to (B.2b), (B.6) for small values of \t\, as we will check below. 

There arc other possibilities for the t dependence of the rho residue. We could have chosen solution 1 
of Rarita et al.;'^^' or we could have taken a Veneziano-type formula with a diffractive factor included, writing 

p{s,t) =13 (i)(s/s)"''^*^; 
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Figure 16. The t-dependence of the different Regge residues, p{t)/f3{0) (refs. 17, 18): Rarita 
et al., solution 1 and solution la: the fits from Rarita et al. (assuming equal slope for ttN and 
tttt). Froggatt & Petersen: from fits to tttt dispersion relations at fixed t. Palou & Yndurain; 
diff. factor: actually, from the Veneziano model, including a diffraction factor e*"', 6 = 2.4 GeV, 
as in Eq. (B.7). 

The residue given in (B.6) would lie between "Palou & Yndurain; diff. factor" and "Rarita et 
al., solution la" (closer to the second). 

Another possibility is afforded by the t dependence obtained by Froggatt and Petersent^^l from an analysis 
of TT+TT" dispersion relations at fixed t. These authors take 

"'("'"'("'Sit' (^^^' 

with t-^ = -2.2 GcV"^ 6fp = 0.8 GeV"^. This gives a slope for the residue of about half the value of 
Solution la of Rarita ct al. 

The values of Pp{t), both /3p(0) and t dependence, are given by Rcggc fits with loss precision than 
other quantities; a few determinations (including the ones mentioned here) are shown in Fig. 16. The 
variation from one estimate to another in the range necessary for Roy equations, —t ~ 0.6 GeV^, is very 
large, and even at [tj - 0.1 GcV^ there are noticeable differences. 

Fortunately, we can at least check the correctness (within errors) of the values of f3p{t) given in ref. 5 
at t = and t — 4Mj ~ 0.08 GeV^, and even improve them. In fact, apart from the sum rule (5.2), the 
fulfillment of the dispersion relation (3.5) fixes with good accuracy /3p(0), and combinations of both values 
/3p(0) and /3p(4M^) enter in the evaluations, using the Froissart-Gribov representation, of two quantities, 
ai, bi, that can also be determined by other methods, independent of the Regge assumptions. The ensuing 
equality of the two estimates each for ai, bi given in Eqs. (6.2) and (6.4) fixes Cp, as we will see below. Thus, 
if, for example, we had taken the Veneziano-type t dependence of Palou and Yndurain^^^l (i.e., as in (B.8) 
without the exponential factor e^*) we would have obtained &i ~ 3 x 10~^ quite incompatible with 

determinations from the pion form factor or the sum rule (6.8). Even if we include a diffraction factor as in 
(B.8), the values of ai, 6i following from the Froissart-Gribov representation vary by more than 1 a. 
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One can in principle fix the parameters for tttt scattering from other processes profiting from factor- 
ization, which allows one to improve the poor precision of the tttt data. This is all right for the P, P' poles 
because here we can use data from ttN and pp+pp, scattering which are very precise and to which only these 
two poles, plus the rho pole for ttN, contribute. The situation is different for the coupling of the rho pole 
to TTTT. Direct fits to tttt data arc not precise and, to use factorization, one has to incorporate (besides ttN) 
proton, antiproton and neutron scattering data. The last presents important systematic errors because the 
data have to be extracted from scattering on deuterons. The large number of trajectories that contribute 
to these processes also make the analysis unreliable; for example, some existing fits have very likely their 
rho trajectory contaminated by pion exchange, which is not negligible at the lower energy range. Because 
of this, we will now give a new determination of /3p(0) and dp independent of Regge fits. 

To obtain pp{0), we rewrite the dispersion relation (3.7) as 

R.f<.=..(.,„,_ Hizi^ep.p. /'"•"•^'d.' /°'fr'''-'°) 

74M2 (S'-S)(S' + S-4M2) ^^^^^ 

2s-4M2 ^ . ImF(-f*=i)(s',0) 

d.s 



'l.42^GeV^ (,s'-5)(,s'+5-4Af2)- 

We evaluate the left hand side for values of s^^^ from threshold to 925 MeV, at intervals of 25 MeV, using 
low energy experimental data, via the parametrizations^ given in Sect. 2. In this region, s^/^ < 0.925 GeV, 
the uncertainties due to the middle energy SO and P waves are not very important. If we write in the right 
hand side of (B.IO) 

ImF(^*=i)(s',0) = /3p(0)(s7s)"''W, s > 1.42^ GeV^ 

and fix ap{0) = 0.52 from factorization or deep inelastic scattering, we may fit (3p{0) and find a very precise 
number, 

/?p(0) = 1.06 ±0.10. (B.ll) 

This is quite compatible with the value /3p(0) = 0.94 ±0.14 that we obtained from fits to tttt cross sections in 
rcf. 5. The result (B.ll) is robust; if we double the errors of the parameters for the P and SO waves between 
1 and 1.42 GeV, where (as discussed in the main text) one has larger uncertainties, and repeat the fit, we 
get the same value. 

As for Cp, we can use the values of ai, bi obtained from the analysis of the pion form factor, 1^1 
together with the Froissart-Gribov representation of the same quantities. We rewrite Eq. (6.1) as 



4M,r(l + 3/2) isM^ 

/^r{2) f°° inijr(/t=i)(g^4M2) 



8M^r(l + 3/2) . 

'1.422GeV2 



/ 

Jl.' 



(B.12a) 



For the 6;, I =odd 

^r{l + l) /-i-^a^GeV^ j 4 Im (s, 4M^) (Z + 1) Im F(g, 4M^) 
' AM^r{l + 3/2)J^M- M (s-4M2)s'+i 



^^/nr{l+l)_ r°° [ 4 Im (g, 4M^) (/ + 1) Im F{s, 4M^) | (B.12b) 

~4M,r(Z + 3/2) A.42.GeV^ '{ {s-4M^)si+^ J' 
ImF(^*=i)'(s,i) = dImF^^'=^\s,t)/dcose. 



® Of course, before improving with dispersion relations, or the reasoning would be somewhat circular. Alternatively, 
we could fit simultaneously the parameters for phase shifts and inelasticities and the Regge parameters of the rho, 
using all three dispersion relations; the method, quite complicated, would not give results substantially different 
from what we get making independent fits. 
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The left hand sides in (B.12) are evaluated using also the experimental phase shifts and inelasticities as fitted 
in Sect. 2 of the present paper and we take, in the right hand side of (B.12a), 

IniF(-f*=i)(s,4M2) = /3p(4M2)(s/s)«p(4M^)^ ^ > ^^ 42^ GeV^ 

Pp{t) ~ 6f,{0)\l + c^jt]. An analogous formula may be written for the right hand side of (B.12b). The 
sum rule for bi gives the more precise constraint for Cp for two reasons: first, it depends directly on the 
derivative of lmF^^*^^^s,t) with respect to t {cos6), hence Cp appears explicitly; and, secondly, the low 
energy contributions to the left hand side in (B.12b) cancel to a large extent, so 61 is given almost exclusively 
by the high energy, Regge formulas; cf. ref. 6. 

We can fit, with these sum rules for ai, 61, plus the dispersion relation for It = 1, the parameters 
/3p(0), Cp simultaneously. As was to be expected, the value of (3p{0) is unchanged and, for Cp, we find the 
values 

_ / (6.0 ± 2.8) GeV"^ from ai 1 on 

''''~\(4.7±0.8)GeV-^ from 61. ^ 

This is compatible with what one has from solution la of Rarita et al.,'-'^^! Cp = (6.6 ± 0.5) GeV~^ in the 
sense that they produce similar values for (3p{4:M^): the decrease in Cp is compensated (within errors) by the 
increase in /3p(0), Eq. (B.l), compared with what we used in ref. 5. Our results for Cp thus justify the use of 
the formulas of Solution la of Rarita ct al.,[^^l for small values of \t\ < 0.1 GcV, as was done in ref. 6, and 
the assumption that (in that range) they are also valid for tttt scattering, within errors. 

We can also include in the fit the information from tttt scattering data (as discussed in ref. 5), since 
it is compatible with what wc found now, and the sum rule (5.2), evaluated with the parameters found in 
Sect. 2 here for the low energy piece. In this case we find what we consider the best results, 

/3p(0) = 1.02 ± 0.11, Cp = (5.4 ± 0.9) GeV"^. (B.14a) 

This implies 

dp = (2.4 ± 0.9) GeV-^ (B.14b) 
Eqs. (B.14) provide us with the central values given in (B.6) for dp, (3p{0). 

Appendix C: On experimental phase shifts in the range 1.4 GeV ~ s^/^ ~ 2 GeV 



As is well known known, as soon as the center of mass kinetic energy in a hadronic reaction, Skin, cm., 
increases beyond 1 GeV, inelastic processes become important and, for i?kin.c.m. ft 1-3 GeV, they dominate 
elastic ones. This is easily understandable in the QCD, ladder version of the Regge picture, as discussed 
in ref. 6; and indeed, it is verified experimentally in the hadronic processes ttTV, KN, and NN where, for 
£^kin,c.m. > 1-3 GcV, the elastic cross section is smaller than the inelastic one and, for -Bkm.c.m. > 1-5 GeV, 
the elastic cross section is less than or about one half of the inelastic one. There is no reason to imagine 
that TTTT scattering would follow a different pattern. In this case (large inelasticity), and again as mentioned 
in refs. 6, 7, it can be proved theoretically!^^! that there is not a unique solution to the phase shift analysis: 
some sets of rjs and Ss may fit the data: but so would others. 

In spite of this, the Cern-Munich experiments^ have produced a set of phase shifts and inelasticities 
which go up to ~ 2 GeV. These have been used as input in certain theoretical analyses, notably in 
those of refs. 1,2. Unfortunately, such phase shifts and inelasticities are very likely to diverge more and more 
from reality as s^^^ becomes larger and larger than 1.3 GeV ~ 2M^ + 1 GeV. This is suggested, besides the 
theoretical reasons just mentioned, by the fact that the Cern-Munich phase shifts and inelasticities contradict 
a number of physical properties they should fulfill: we will here mention a few. 

As we have remarked above, one would expect dominant inelastic cross sections above s^^'^ ~ 
1.5 GeV; but the Cern-Munich elastic cross sections are larger or comparable to the inelastic ones up 
to sV2 := 2 GeV. What is worse, the tt+tt inelastic cross section of Cern-Munich, alone of all hadronic 
cross sections, decreases as the kinetic energy increases between 1 GeV and 1.7 GeV, as shown for example in 

Hyams et al. and Grayer et ai., refs. lla,b. 
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Fig. 7 in the paper of Hyams et alJ^^"! As happens for jN, -kN, KN or NN scattering, and indeed as is also 
seen in other experimental analyses for tttt scattering, I**! wc expect a levelling oif of the total cross section for 
-E^kin.c.m. > 1-3 GeV: but the Cern-Munich total cross section for tt+tt" scattering decreases roughly like 1/s 
up to 2 GeV. In fact, the Cern-Munich elastic cross section agrees well with that found in other experiments 
(such as those in ref. 4), but their total cross section is smaller by a factor ^ 2 than those in refs. 4, at 
the higher part of the energy range: the Cern-Munich results certainly misrepresent the totai cross section. 
Prom all this it follows that the Cern-Munich phases and inelasticities are not reliable in that energy region. 

Secondly, the combination of S and rj for both P and SO waves at energy > 1.8 GeV is incompatible 
with QCD results for the electromagnetic form factor of the pion, and also with Regge behavior, which 
requires all phases to go to a multiple of tt at high energy, and the real part of the partial wave amplitudes 
for isospin 0, 1, to be positive. The phase Si is also incompatible with QCD results for the electromagnetic 
form factor of the pion. In fact, the phase of this form factor behaves like 

^Form factor(s) TT ( 1 + 5— ^ ) , S > yl^; t 

\ logs/tj 

A is the QCD parameter (see, e.g., ref. 20). One may take this to hold for s ^ 3 GeV^ (s^/^ > 1.6 GeV). 
If one had negligible inelasticity for the P wave somewhere in the region 1.6 GeV < sV2 < 2 GeV, as the 
Cern-Munich data seem to imply, the form factor and partial wave would have the same phase at such 
energies, and thus the same behaviour should hold for 6i{s), c^Porm factor (s). But the phase which the Cern- 
Munich experiment gives clearly contradict this behaviour around 1.8 GeV; there the Cern-Munich phase 
5i stays consistently beJow tt, while, as just shown, it should be above. 

Thirdly, for the process tt^tt^ — > tt^tt^, we would expect large inelasticity as soon as the production 
TT^TT^ p^p^ is energetically possible, and therefore large inelasticity for the S2, D2 waves for s^/^ » 2Mp. 
This, in particular, occurs in any theoretical model. It is therefore not possible to extrapolate these phases 
above ~ 1.45 GeVas being elastic. In fact, the extrapolation used in refs. 1, 2 for the D2 wave is clearly 
incorrect above -^1.4 GeV as the corresponding |(52^''(s)| grows linearly with s while, from Regge theory 
(and also from low energy fits, see Fig. 7), one expects it to go to zero. Thus, besides the general problem 
for the cross sections we have individual problems for each of the SO, S2, P and D2 phases. 

Finally we would like to mention that both the Regge picture and the values of the experimental 
cross sections for all hadronic processes indicate that the number of waves that contribute effectively to the 
imaginary part (say) of the scattering amplitudes grows with the kinetic energy as Ey^I A, for ii^kin upwards 
of 1 GeV. Wc thus expect 2 to 3 waves (for fixed isospin) at £^kin ~ 1 GeV, and almost double this, 4 or 5 
waves, at ii^kin ~ 1-7 GeV. In fact, for tttt scattering at this energy, the contribution to the total cross section 
of the F wave is as large as that of the P wave, the DO wave is as large as the SO wave and the contribution 
of the D2 wave is actually larger than that of the S2 wave. The partial wave series with only two waves per 
isospin channel is not convergent, and the approximations, like those of refs. 1, 2, 11, that neglect higher 
waves at such energies have yet another reason for being irrealistic. 

Because of all this, it follows that use of the Cern-Munich phases and inelasticities must lead to 
a rather distorted imaginary part of the tttt scattering amplitude above s^l"^ ~ 1.4 GeV. It is thus not 
surprising that authors like those in refs. 1, 2, 15, who fix their Regge parameters by balancing them above 
2 GeV with Cern-Munich phase shift analyses below 2 GeV, get incorrect values for the first. ^ 



We would like to emphasize that what has been said should not be taken as implying criticism of the Cern-Munich 
experiment which, for s^^^ ;$ 1.4 GeV, produced what axe probably the best determinations of phase shifts and 
inelasticities. Above 1.4 GeV, they did what they could: it is for theorists to realize that this was not enough to 
produce acceptable phase shifts and inelasticities at these higher energies. 
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